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ABSTRACT 


In the summer of 1969 a line of seven magnetic 
stations was established in western Canada between the 
geomagnetic latitudes of 58.5° and 77.0°N and within 2° 
of 302°E corrected geomagnetic longitude. This thesis 
gives an analysis of the general, spectral characteristics 
(including the three-dimensional polarizations and station- 
to-station phases) of Pc 4 and Pc 5 geomagnetic micro- 
pulSations recorded at these stations in 1969 and 1970. 

A change in the horizontal sense of polariz- 
ti ON Ot, Lie MICroDUlsSaLicns , occurring around rlUU? to 
1200 local geomagnetic time, indicates that Pc 4's and 
Pe 5's are generated through the development of Kelvin- 
Helmholtz instabilities at the magnetopause. The power 
Spectra and polarizations Or the InpcropulsatiTons @xnibi t 
marked latitudinal changes, implying that much of the 
micropulsation energy is distributed in resonant oscil- 
lations of magnetic shells. The existence of distinct 
Pe 4 and Pc 5 spectral bands may be the result of prefer- 
ential coupling of energy into magnetic shells near the 
auroral oval and near the plasmapause, with Pc 5's 
CCCURYING Drifiarl ly an tne auroral regions and Pc 44s 


near the plasmapause. 
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1. INTRODUCTION 


Tole Geomagqneticimicropulsations, deninition. and classi — 


fication 


Geomagnetic micropulsations are transient, ultra 
low frequency oscillations of the earth's magnetic field. 


3 Hz) often have 


Low frequency micropulsations (*10° 
amplitudes of hundreds of gammas at the earth's surface 
although high frequency micropulsations (*1 Hz) seldom 
have amplitudes larger than a fraction of a gamma (y). 

The structure of the pulsations may be very irregular and 
ragged or almost sinusoidal (Figure 1). In some cases the 
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iS continue over 
pulsation trains lasting for hours. At other times the 
pulsations disappear after only one or two oscillations. 
It is now commonly believed that continuous micropulsations 
originate within the magnetosphere as hydromagnetic waves, 
although there is considerable uncertainty about the energy 
sources and modes of propagation of these waves. 
Micropulsations are normally classified in three 
characteristic ways, by their morphological properties 
(frequency, wave form, etc.), by their correlative properties, 


and by their genetical properties [Jacobs,1970]. At present 
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Fig. 1. An example of low frequency, continuous micro- 
pulsations. This event was recorded in the interval 
1430-1600 UT, 1970. 
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both the correlative and genetical properties form a poor 
basis for classification in that correlative properties 

have been observed in only a few cases and, as mentioned 
above, the specific origins of many types of micropulsations 
are uncertain. A classification system based on the 
frequency of the oscillations and the wave form of the 
pulsation train is now commonly accepted. In this system 
regular, continuous pulsations are denoted by the term Pc 
and irregular, impulsive pulsations by the term Pi. Both 
the Pc and Pi classes are divided into frequency ranges 
based on other physical and morphological properties [Jacobs 
Arar Uwe Loo. ehdaviles leGLVEGS cawLiS LRORMEneSeRClaSSini Cations. 
Hultqvist [1966] gives a more comprehensive description of 


the different classification systems. 
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The experiment outlined in this thesis was designed 
to study the morphological features of continuous micro- 
pulsations in the spectral range 1-30 mHz (Pc4 and Pc5) 
at a series of ground stations set at different latitudes 
along a geomagnetic meridian. For comparative purposes the 
analysis includes a number of nighttime Pi micropulsation 
events. We have reduced the data by means of a generalized 
SDECULdMEanalLySismeOf Ener MiCropuisdatnon  silgndis-= in tins 
method of analysis the autospectra and cross spectra are 
arranged in the form of 3x3 matrices, based on the three 
spatial components, and from these matrices the intensities, 
polarizations, and station-to-station phase characteristics 
of the micropulsations are calculated. This type of 
analysis allows a more objective and complete approach and 
puts the results in a format which is more compatible with 
the literature on the theory of micropulsations. 

The analysis indicates that much of the energy of 
dayside Pc 4's and Pc 5's, and perhaps some of the energy 
of nightside Pc's,comes from Kelvin-Helmholtz instabilities 
at the dawn and dusk magnetopause. At least some of tnis 
energy is coupled into resonant oscillations of magnetic 
shells within the inner "dipole" region of the earth's 
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spectral bands may be the result of preferential coupling 
of energy into magnetic shells near the auroral oval and 
near thempiliasmapausers with Pe 5 'ssocctinning primarily. in 


the auroral regions and Pc 4's near the plasmapause. 


1.2 Review of the observational results 


Sincewoteuart | | Coll enirstereportedemi cropulsations:. 
a great diversity of papers has appeared in the liter- 
ature. The largest increase in the output of papers came 
aDOUCe eet iewmlates O50 ss at teretnerlaGwyen sinethe [9602s 
satellite data first appeared but, as yet, the number of 
reports of micropulsations is limited. The purpose of this 
review is to discuss very briefly the observations which 
are most pertinent to the data presented in this thesis. 
The first section deals with ground-based observations of 
Pc 4's and Pc 5's and the second section with recent satel- 
ltemobsenvations . srotteccyva [196/ ln Gampbel? \[196/1, 
Satto [1969], Dungey and Southwood [1970], and Jacobs 


[1970] have given more comprehensive reviews in these areas. 
(A) Ground-Based Observations 


Both Pc 4's and Pc 5's have maximum amplitudes 
and occurrence frequencies in the auroral zones (6525708 
geomagnetic latitude) [Jacobs and Sinno, 1960; Obertz and 


Raspopov, 1968; and Hirasawa, 1970] although Kato and Saito 
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[1962] have noticed a second Pc 4 peak at 50°. The amp li- 
tudes of Pc 5's often reach hundreds of gammas in the auroral 
region, decreasing to about 1/5 of this value at 80° and 
50°. Pc 4 amplitudes, on the other hand, are seldom greater 
than l10y. The longitudinal extent of both Pc 4's and Pc 5's 
far exceeds the latitudinal extent and follows roughly an 
east-west direction along the auroral zone [Eleman, 1966]. 
Obertz and Raspopov [1968] observed that the characteristic 
Wonda Gudin alee Xten Grote Gamo. Suis 60° with most activity 
centred between 0600 and 0900 LMT. 

The diurnal variation of the occurrence frequency 
Ofeauroral zone; Pe 5°s dififiersstrom that of. midlatitude 
Pica Seis (40°-55%). In the auroral zones the activity peak is 
sharply defined and centred near 0400-1000 LMT, with a 
small second peak in the afternoon [Saito, 1964; Hirasawa, 
1970] (see Figure 2). In sub auroral regions the occurrence 
frequency is double peaked with a second peak in the after- 
noon [0Z' 1963; Saito, 1964]. In midlatitude regions Pc 5's 
occur with equal probability over most of the day. 

The amount of Pc 4 activity usually has a broad 
maximum in the daytime at midlatitudes [Satto, 1964] whereas 
the maximum is at 0300-0600 LMT in the auroral zones [Saito, 


1964; Stuart and Usher, 1966; Annexstad and Wilson, 1968]. 
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from a single maximum at high latitudes to a double 
maximum at low latitudes (after Saito [1964]). 
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In general both Pc 4's and Pc 5's are most sinusoidal in 

the morning and early afternoon, the greater part of the 

evening and nighttime activity having only a few damped 

cycles in the wave train [Jacobs and Wright, 1965]. 
Statistical studies by Obayashi and Jacobs [1958], 

O7utz se (1963) and@o7* 11963] “showethat the periods of Pc 4's 

and Pc 5's increase with increasing latitude (Figure 3). 

O1' concluded that the average period T is given by the 


equation 


T= C-sec“é 


where 6% is geomagnetic latitude and C is a constant 
depending on diurnal and solar cycle effects. Normally, 
however, a given micropulsation event has the same frequency 
at all latitudes [#ZZis, 1960; Obertz and Raspopov, 1968]. 
The latitude-dependent Pc's noted by Voelker [1968] are 
exceptions to this rule. These micropulsations exhibit a 
systematic increase in the dominant period of the H component 
with increasing latitude. This particular type of event will 
be discussed in more detail in the analysis of data in this 
thesis. 

Saito [1969] postulated that although the frequency 
of a specific micropulsation event is the same at all latitudes 


(Sea profiles. da sande bo in Figure 4), the statistically 
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Fig. 3. The latitudinal dependence of the periods of 
P Coteasmo oli CusrCeeones ue WeLAV a [UCSHAR/i ea rem tnemequia tor iva, | 
radii of the field slines in units! of earth radii (after 


Saito [1969]). 
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Fig. 4. Latitudinal and period dependence of Pc 4 
and Pc 5 amplitudes. The curves a@ and b_ represent 
the latitudinal dependence of the amplitudes of Pc 4's 
and Pc 5's respectively. When the geomagnetic activity 
increases, the curves a and b move along line e 
COMM OWermeld titudest ee Peaks ofA a BigiandeGibe long) toyP ca, 
Pc 5, and giant pulsations respectively (after Saito 
[1969]). 
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determined frequencies can decrease with increasing lat- 
itude. He suggested that if the intensity peak of a Pc 4 
or Pc 5 micropulsation train moves southward, due to 
increasing magnetic activity, then the frequencies of the 
pulsations increase (i.e. the peaks A and B move along 
line c in Figure 4). The analysis of Obertz and Ras- 
popov [1968] supports Saito's model. 

In general both Pc 4's and Pc 5's have higher 
frequencies during the day than at night. At midlatitude 
Stations the Pc 4 frequency reaches a broad maximum centred 
near 1200 LMT [Christoffel and Linford, 1966; Saito, 1964]. 
The maximum frequency, occurring during the day, is 30-50 
mHz and the minimum frequency, occurring at night, is 10-15 
mHz. Htrasawa and Nagata [1966] and Nagata and Fukunishi 
[1968] found that the midlatitude Pc 4's actually occupy 
two distinct spectral bands, one centred at 15 mHz and the 
other at 30 mHz. The 30 mHz band is most evident during 
the day. Both bands show pronounced diurnal variations in 
frequency (Figure 5). During periods of low magnetic activity 
(Kp 0v1) the highest frequencies occur around 1500 LMT and 
during periods of high activity (K, 5v7) around 0700 LMT. 


In the auroral zones, Pc 5's have their highest frequencies 


around 0600 LMT [Satto, 1964; Hirasawa, 1970] but in sub 
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Fig. 5. Average sonograms of geomagnetic micro- 
pulsations recorded at Kakioka. 


(a) Quiet conditions. Ky = (bh ay eR 
(b) Slightly disturbed conditions, Ky 
(after Nagata and Fukunishi [1968]). 
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auroral and midlatitude regions they have their highest 
frequencies around 1200 LMT [xttamura, 19633 Kato and 
Suvto ew t90c |. 

Both Pc 4's and Pc 5's show polarization char- 
acteristics which change with latitude and time. Normally 
the vertical component is the smallest but it does reach 
appreciable amplitudes in the auroral regions [EZeman, 1966]. 
A complete analysis of the polarization characteristics of 
micropulsations requires the determination of the projec- 
tions of the polarization ellipse on three planes. Unfor- 
tunately, this has only been done for a few events [Paulson 
et al., 19653; Wilson, 1966; Annexstad and ele? 1968]. 

To date most analyses have dealt only with the horizontal 
projections of the polarization ellipses. 

In northern high latitude regions (> 60°), the 
NoOrizonca | components O01 Pores S=arempolarized in a CC 
(counterclockwise)sense in the morning (looking down on 
the earth), with 8-10 hours of predominantly CW (clockwise) 
polarization during the afternoon and evening. The transition 
between the morning and afternoon sectors typically occurs 
around 1100 local geomagnetic time [Kato and Utsumt, 1964; 
see Figure 6]. The morning section of CC polarization is 
the clearest and the afternoon section most confused, having 


an appreciable amount of both CW and CC polarization. These 
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Point Barrow July 1957-June 1958 
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field vector in the horizontal plane (after Kato and 


Utsumt [1964]). 
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diurnal polarization characteristics are reversed in the 
southern auroral regions [Satto, 1969]. In middle latitudes, 
the diurnal polarization characteristics of Pc 5's are more 
complex, with four sectors alternating from CW to CC [Satto, 
1964]. 

The diurnal polarization pattern of high latitude 
Pc 4's matches the Pc 5 pattern (see Figure 7). One strik- 
ing feature of Figure 7 is that the sense of polarization 
Of seCuz.. sus) (30-100MmHz) i1S7alWayomonposite to that of 
PCO seeeCeeasmande PC.o4S NavemGnmpoldrization in. tie 
morning with a pronounced change to CC polarization in the 
afternoon. In northern midlatitudes Mather et al. [1964] 
found CW polarization for daytime Pc 4's and CC polarization 
for nighttime Pc 4's. This observation agrees with the 
southern, midlatitude Pc 4 polarizations determined by 
Christoffel and Linford [1966] (i.e. CC during the day and 
CW during the night, looking down at the earth). 

In addition to diurnal changes in polarization, 
Pe 455 "s exhibit’ pronounced latitudinal changes in their 
sense of polarization. Kaneda et al. [1964] and Obertz 
and Raspopov [1968] reported that Pc 5's at Barrow (68.6°) 
often had CW polarizations in the horizontal plane when 


the polarizations at Sitka (60.0°) and College (64.7°) 
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Figs 7. Average senses of polarization in the 
horizontal plane of micropulsations in the spectral 
range 3.3 - 100 mHz. The shaded regions marked C 
have predominantly clockwise polarization and those 
marked CC have predominantly counterclockwise polar- 
ization. The unshaded regions have no predominant 
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sense of polarization (after Rankin and Kurtz [1970]). 
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were CC. An appreciable latitudinal phase shift in the H 
component accompanied the polarization reversal. In 
addition, the opposite sense of polarization appeared only 
when the centre of activity was between College and Barrow. 
Eleman [1966] found similar latitudinal changes in the senses 
of polarization of Pc 3, 4's, and Sakurai [1970] found 
Si ntinler chia Wace SeCSatf Oe Pagiceuse 

Conjugate observations of pulsation waveforms and 
polarizations are important since they directly relate 
micropulsations to the magnetosphere. Studies by Nagata 
et al. [1963] and Jacobs and Wright [1965] showed similar 
Pe 5 pulsations occurring simultaneously at conjugate points. 
Theyanound.. ini taddi tion .ethat: thebsensesaof polarization 
were conserved along the geomagnetic lines of force. Annex- 
stad and Wilson [1968] observed the same characteristics 


for higher frequency (10-15 mHz) Pc 4's. 


(B) Satellite Observations in the Magnetosphere 


To date few observations of micropulsations have 
been made by satellites and even fewer observations have been 
related to ground-based data. Pioneers I and V and Explorer 
6 recorded some of the earlier observations of pulsations 
in the magnetosphere [Coleman et al., 1960; Sonett et al., 


1962; Judge and Coleman, 1962]. Patel and Cahill [1964] 
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and PateZ [1965] reported irregular 5.5-8.3 mHz (120- 

180 sec period) transverse oscillations recorded simul- 
taneously both in tne magnetosphere by Explorer 12 and 

at the ground on the same magnetic field line at College, 
Alaska. The waves needed 1 1/2 min to travel from the 
Sateiiite sto the oqround, 9Pategmand Cant ll) also found 

that the polarizations, frequencies, and amplitudes of 

the pulsations at College agreed with those at the satel- 
lite. In addition Patel found a local-time-dependence for 
the polarization at the satellite which was substantially 
in agreement with ground-based Pc 5 observations (i.e. CC 
before 1100 local geomagnetic time, CW after). Dungey and 
Southwood [1970], in an analysis of magnetic field data 
from Explorer 33, determined the same diurnal polarization 
rule for fluctuations inside the magnetosphere near the 
magnetopause (Figure 8). 

One of the first reports on observations of con- 
tinuous pulsations in the magnetosphere was made by Cummings 
et al. [1969] using records from the geostationary satellite 
ATS 1. These pulsations occurred during magnetically quiet 
periods, mostly in the interval 0400-1600 LMT. The oscil- 
lations were very sinusoidal and many continued for 10-400 


min before disappearing (Figure 9). The observed frequencies 
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EARTHSIDE 
POLARISATION 


o CLOCKWISE 
o ANTICLOCKWISE 


Y (Re) 


Eldeo, weOrbitsS=Of -EXDINGeruooeDrojected on the 
ecliptic plane. The magnetopause crossings and the 
polarizations of the micropulsations are marked by 


open or closed circles (after Dungey and Southwood; 
C1970} )e 
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Fig. 9. An example of the transverse oscillations 
observed eat Atsi lo The starting time (UT) s in the 
upper right corner of the diagram (month, day, hour, 
year). The D component is positive eastward, the 
V component is positive outward from the earth, and 
the H component is positive northward (perpendicular 
to the equatorial plane) (after Cummings et al. [1969]). 
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occupied two broad bands, one centred at 10 mHz (102 sec 
period) and the other at 5.3 mHz (209 sec period). Most 
of the pulsations were elliptically polarized in a plane 
transverse to the main magnetic field. More recently 
Barfield et al. [1971] have reported an observation of 
compressional Pc 4 micropulsations at ATS 1. Both the 
magnetic field intensity and energetic electron fluxes 
exhibited sinusoidal modulations at 9.4 mHz. These oscil- 
lations lasted for nearly seven hours. The largest 
pulsations reach 8y (peak to peak) in the direction of 

the main magnetic field, with almost no magnetic variations 


transverse to the main field. 


InSeeEXCT Cation mecnanismiss proposedaror continuous micro-= 
pulsations 
Continuous micropulsations are now thought to 
originate as hydromagnetic (hm) waves in the magnetosphere 
but there is still considerable uncertainty concerning the 
sources of energy which excite these waves. Some of the 


current theories are those based on 


(a) gradients in particle and ionization 
densities [Swift, 1967]. 
(b) non neutral plasma distributions [Wishida, 


1964; Kimura and Matsumoto, 1968]. 
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(c) resonant interactions of particles and 
waves [Parker, 19613; Dungey, 1965; South- 
wood et al., 1969; Dungey and Southwood, 
1970] 

(d) surface wave motions at the boundaries 
between two regions of different plasma 
concentrations [Dungey, 1954; Sen, 1965; 
Atkinson and Watanabe, 1966; Southwood, 


1968]. 


These mechanisms are discussed in some detail in the follow- 
ing sections. 

The energy in a gradient instability is derived 
from the free energy of the inhomogeneous plasma [Swift, 
1967]. The instability results in the generation of electro- 
Static waves in the magnetosphere which couple to electro- 
magnetic waves in the ionosphere. The electrostatic waves 
have CW polarization (looking in the direction of the 
magnetic field) in a region of enhanced plasma density and 
CC polarization in a depleted region. 

Instabilities due to non neutral plasma distrib- 
utions are excited by an electron or ion beam (such as 
precipitating auroral particles) passing through some part 
of the magnetosphere [Kimura and Matsumoto, 1967]. The 


frequency of this instability depends only on the magnitude 
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of, the’ excess charge and. can be as, low as» the Pc 5 spectral 
band. A negative excess charge causes an instability in 
the Alfvén mode whereas a positive excess causes an instab- 
ility in the modified Alfvén mode. 

The most natural wave-particle instability for 
Pc 4's and Pc 5's is one involving bounce resonance [Dungey 
and Southwood, 1970]. Southwood et al. [1969] have shown 
that this resonance can occur in the magnetosphere under 
certain, very specific conditions. They found that quasi- 
transverse hm modes of large m _ can be excited by the 
exchange of energy with energetic particles. (The azimuthal 
variations of the wave are taken as e Md where » is the 
azimuthal angle and m is an integer.) The change in the 
energy of a particle is directly proportional to the mag- 
netic shell coordinate (L) of the particle and is independent 
of the particle's energy or pitch angle. The dominant para- 
meter in energy exchange is the tilting of field lines and 
under certain conditions the growth rate of the waves can 
be very large. 

The surface wave or Kelvin-Helmholtz instability at 
the magnetospheric boundary has long been a popular mechanism 
for the excitation of low frequency micropulsations [Dungey, 


1954]. The surface waves are generated near the sun side 


of the equatorial magnetopause by the tangential flow of 
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the solar wind around the magnetosphere, flowing to the west 
on the morning side of 1100 local geomagnetic time and to 
the east on the evening side [Atkinson and Watanabe, 1966]. 
This tangential flow of the wind forces fluid elements in 
the magnetosphere to rotate in the manner shown in Figure 10. 
Since the magnetospheric plasma in this region is a good 
conductor the magnetic field is "frozen" into the plasma 

and is forced to rotate with it. Thus we expect the polar- 
PZautOnmcONnttagunatuonmOt Figure sluilnewiucn CG poOlarizatyi on 
before 1100 local geomagnetic time and CW polarization after 
1100. .These polarizations agree with both satellite and 
ground-based observations [Dungey and Southwood,1970]. 
Southwood [1968] deduced that in fact the magnetospheric 
boundary can be unstable. He showed that in equatorial 
regions, away from the subsolar point, the first growing 

hm modes propagate perpendicular to the geomagnetic field 
and are circularly polarized in a plane perpendicular to 


the main magnetic field. 


1.4 Normal modes of oscillation in the magnetosphere 


The theoretical treatment of hm modes in the 
earth's magnetosphere is extremely difficult. The problem 
is obscured by mathematical difficulties, by the complexity 


of the resonant cavities within the magnetosphere, and by 
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Fig. 10.. Motion, ofsfluid elements in a free.surface 
wave (after Atkinson and Watanabe [1966]). 
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MORNING AFTERNOON <! 


Fig. ©! LAgeliie senses of polarization@oT geomagnetic 
micropulsations induced by a Kelvin-Helmholtz instab- 
ility at the magnetopause. 


a _ 0 
i 
Daan F 7 = i 
7 - 7 
' 
Ss, in 
3 
2 
4 
7] 4 ” 
i 
SN MAL, > th ew © wat. Saf w oa=e oe yak wie te > - eka SGA aie ap eee wench | Ctra 
A 
y } 
fj ‘, . 
’ 
§ 
) e 
‘ : 7 
- > f 
poe j 
: j 
f 
. f 
Ce™ St 
i we” t 7 ros 
_ 
7 
: i o 
512 eppedo: Ai Jena oe 1v ae SP i 2. ay 4] 


; -UGvtin: anion 1 @r- iV vans # ue “ol PUGS TT. F I 
ae | Hu AGoxelnsin “ont 


: Raeic ) 
1 o< vs Dive 


26 


uncertainties in the correct boundary conditions. Up to 
the present, most theoretical treatments have considered 
the resonant modes of an axially symmetric magnetosphere 
[Dungey, 1954; Jacobs and Westphal, 1964]. If we further 
assume that the disturbance field is axisymmetric, then 

two simple hm modes exist. The toroidal mode has the 
perturbation magnetic field polarized in an azimuthal 
(east-west) direction and corresponds to a torsional oscil- 
lation of a magnetic shell. The second mode, the poloidal 
mode, corresponds to a compressional oscillation of the 


whole magnetosphere. Normally, however, these modes are 


not good approximations to actual conditions sinca most 
observations show that the longitudinal (azimuthal) extent 
OfaPords Ss ismlimited. Mavncen [1954051967] and fadosk 
[1967] have discussed very briefly some of the characteristics 
of highly asymmetric perturbations. 

The basic configuration of the magnetosphere and 
1ts resonant cavities is depicted in Figure 12. The earth- 
ionosphere boundary is usually considered a perfect reflector 
for low frequency hm waves because of the high apparent 
conductivity of the earth at these frequencies .Dungey [1970] 


has indicated, however, that this assumption is open to some 


doubt. 
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A problem related to the reflection of hm and 
electromagnetic waveswithin the earth and ionosphere is 
the determination of ground level magnetic fields which 
result when hm waves interact with the ionosphere. This 
problem is particularly important if we wish to relate 
ground-based observations to magnetospheric phenomena. 
Consequently the transmission and reflection of hm waves 
in the ionosphere has been considered in detail in a 
number of theoretical discussions, including those by 
Prince and Bostick [1964], Fteld and Greifinger [1965], 
Gretfinger and Greifinger [1965] and Inoue and Schaeffer 
[1970]. 

The plasmapause (Figure 12) is a field-aligned 
density discontinuity which separates the plasmasphere, a 
region of relatively high plasma density (» 100 electrons/cm>), 
from the plasma trough (density » 1 electron/cm°). In the 
plasma trough the number density of electrons decreases as 
rt whereas inside the plasmasphere the density is determined 
by diffusive equilibrium [Angerami and Carpenter, 1966]. 
During periods of moderate magnetic activity (K,, = 2-4), 
the geocentric distance to the dayside plasmapause is about 
3.5 Rp (earth radii) whereas the distance to the nightside 


plasmapause is as great at 5.5 R-. Carpenter [1970] and 
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Chappell et al. [1970] have given detailed descriptions 
of the dynamic behavior of the plasmasphere and plasma- 
pause. 

Near the magnetic poles the geomagnetic field 
loses its approximate dipole shape and the field lines are 
swept back into the magnetotail configuration. This region 
also includes part of the plasma sheet and the electron 
densities are higher than those in the plasma trough 
(0.3-30 Acute Get yeas 1968a))eee ue outer 1imits of 
the geomagnetic field are marked by the magnetopause. In 
the solar direction the geocentric distance to the magneto- 
pause is typically 10 Re although distances ranging from 
have been observed. 
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Most theoretical treatments of magnetospheric 


B 


oscillations use hm equations to describe the perturbation 
magnetic and electric fields. For simplicity and convenience 
we will follow this precedent and use these equations in 

COPSMCIs CUSSION  Sevenethouqghetteisenot, always¥clear that 

the assumptions applicable to hm (or magneto-hydrodynamic) 
equations are justified within the magnetosphere. dHultqvist 
[1966], Holt and Haskell [1965], and Rosst and Olbert [1970] 
have given comprehensive accounts of the conditions necessary 


for the application of the hm equations. 
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The linearized, hm wave equation for a plasma with 


infinite conductivity is 


=> 
> = Va x Va XE Vex Voxee (1) 


(gaussian units) 


where E is the DELLUPD atl OnNwelec trl CamiaLen a. 

Va = B(4mo)7? and Bo is the unperturbed magnetic field. 
The component equations are best analysed by first expressing 
them in orthogonal curvilinear coordinates (xj; Xo» X) with 
the unit vector % in the direction of ee Displacements 


along the three coordinates are given by the relation 


dS = h, dx. (2) 


where the scale factors hs are in general functions of 


X12 Xo, and Xa. In the analysis presented here we assume 
that the system is axisymmetric with respect to x3. (In 
the earth's field, x3, is geomagnetic longitude). If we 
further assume a time-dependence of the form @ Ot then 


the component equations for an axisymmetric system are 


- — 
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Equations (3) and (4) are coupled and consequently are 
very difficult to solve except in a few simple cases. 
Since the system is axisymmetric with respect 


to X93 E. can be represented by a summation of modes. 


oy m imx 
E. = , N. e 3 (5) 


In solving for the mode in which m = 0 (the mode in which 


More precisely, 


the perturbation fields are also axisymmetric), two 


decoupled equations result. 


h 2 
i ale sEolhg%H) ee oe (6) 


2 on (73) 


valwuer eve (9) baw (C)) smPeeaps 


‘ius hres * 3¥ 
s 7 


Be 


Equation (6) describes the toroidal mode of oscillation 
and equation (7) the poloidal mode. The toroidal mode 
appears to be a "guided" mode since there are derivatives 
only along a field line. Equation (6) can be solved by 
using the W.K.B. approximation where applicable [Radoski, 
1966]. By using the W.K.B. approximation and assuming 
that the earth is a perfect reflector, we find that the 


resonant frequencies of this mode are given by the relation 


25 dS, 
of —— = nt test. ey eS igeeee 1D) 
V 
a 
24 


where o, and are the points where the field line 
intercepts the/earth. From this equation it appears that 
w is a function of position and this violates the assumptions 


used to obtain equations (3) and (4). (i.e. we assumed 


that =2- (w) = uwge- .) aruSka [1968] and Kahalas [1969] 
Xs 9X. 
have shown that this dilemma can be resolved by a closer 


inspection of the physical parameters in the problem. Since 
the integral in equation (8) normally increases with in- 
creasing geomagnetic latitude, the expected frequency of 
the toroidal mode must decrease with increasing latitude 


[Dungey, 1954; Jacobs and Westphal, 1964]. 
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The perturbation magnetic field B is found 


using Faraday's Law 


Vane a + a (9) 
More explicitty 
(a, te 88) 
(10) 
cliw) if = 5% ( 33) 3a (P33) oa oe (hoN) 


Equation (10) shows that the magnetic field of the toroidal 
mode is poiarized in tne X 3 direction (east-west) and 
that the poloidal mode has components along the field and 


in the direction. The poloidal mode will not be 


a2 
discussed in any more detail here since this mode has very 
complex analytic solutions. Kato and Akasofu [1955], 
Jacobs and Westphal [1964], and Carovillano et al. [1966] 
have given comprehensive discussions of this mode. 

Another simple but interesting case occurs when 
the disturbance has an extreme longitudinal asymmetry (modes 
of large m). The magnetic disturbance then becomes trans- 
versely polarized in the Xo direction and obeys an equation 


very similar to the toroidal wave equation. Radoski [1967] 


called this mode a guided poloidal mode. 
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Unfortunately the picture presented above is 
drastically simplified. The magnetosphere is clearly 
not symmetric and the approximations used to obtain the 
hm sues BR ai Balas To are open to poneie rable 
debate. A more realistic approach to the problem would 
be to consider the effects of finite conductivity and 
finite Hall currents in the Ohm's law relation and 
possible viscous effects in the equation of motion. 
Inclusion of any of these terms couples the wave equations. 
Even in the simple axisymmetric case, the toroidal and 
poloidal modes are coupled. Magnetospheric oscillations 
with mode coupling have been considered briefly by Hruska 
[1968, 1969] and Kahalas [1969] and in more detail by 


MeClay [1970]. 
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In the summer of 1969 a line of 7 magnetic 
stations was established in Alberta and the Northwest 
Territories of Canada (see Figure 13) to record polar 
substorms and low frequency (1-20 mHz) micropulsation 
activity. In 1970 the station at Penhold was moved to 
Leduc and an additional station was established at Fort 
Reliance to give a total of eight stations. The stations 
lay within 2° of 302° corrected geomagnetic longitude 
[Hakura, 1965]; the southernmost station was at Calgary 


SOc 
ff 


ne et) ; P A ; ms A Fe ae = PS 
(58.5 N) and the northernmost at Cambridge Bay (77 N) 


(see Table 2). Each station used three-component fluxgate 


magnetometers oriented in (H, D, Z) coordinates (Figure 14). 


H (MAGNETIC NORTH) 


D (MAGNETIC EAST) 


Z (DOWN) 
Fig. 14. The magnetic coordinates 
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Data were recorded digitally on 7-track, magnetic tapes 
at a rate of 1 sample/component/sec in 1969, and at a 
rate of 1 sample/component/2 sec in 1970, To time the 
records, each system recorded WWVB directly on tape every 
LOmiina2o mine in O69 andseveryesenresce minein 1970. 
Appendix A6 gives a more complete description of the 
instrumentation. 

The most common source of lost data originated 
from failures in the digital systems. Any large transients 
in line power caused lost data blocks with subsequent data 
gaps and inaccurate timing. This problem was especially 
severe in 1969 since the digital system had to be restarted 
manually after every power failure. An automatic restart 
mechanism was incorporated into the system in 1970. At RELI, 
nearby radio transmitters caused large amounts of digital 
and magnetometer noise. The number of transients in the 
magnetic data from RELI is so great that no data from this 
station are included in the spectral analysis. Other 
difficulties occurring in the magnetometers resulted in 
lost data at SMIT, MCMU, and MENK. At MENK a malfunction 
in the D component amplifier during most of the summer of 


1970 caused sporadic changes in the gain of the channel. 
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TABLE 2 The Magnetic Stations 


Corrected Time Interval 
Station Geomagnetic in Operation 
Coordinates . 
1969 1970 
Aug-Nov May-Nov 
Aug-Nov 


June-Nov 


Sep-Nov May-Nov 


Fort McMurray : : Aug-Nov May-Nov 


Fort Chinewvan : 03. Aug-Novy May-Noy 
Fort Smith : : Aug-Nov May-Nov 
Fort Reliance ~4 : July-Oct 


Cambridge Bay : ’ Aug-Nov June-Oct 


ewe Geoelectric characteristics of the earth in central Canada: 


Associated induction effects 


The geoelectric properties of the earth in central 
Canada are summarized in Figure 15. The electrical properties 
of the substrata in southern Alberta and southern British 
Columbia are well known since many magnetotelluric and geo- 


magnetic depth sounding surveys have been made in these regions 
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Figure 15. Geological and geoelectric regions 
of central Canada 


The numbers in parenthesis on the map are the surface 


resistivities of the regions in units of ohm meters. 


Legend 


E 


The Cordilleran Region, Paleozoic sediments, Tertiary 
volcanics, Mesozoic and Cenozoic intrusives. 


*Surface resistivity 100-1000 Q M. 


Lower and Upper Cretaceous marine and non marine sedi- 
mentary rocks. 


eresistivity 5-400 @ M. 


PROtewOZO1 CaniitUCUSTONS, acCncaic rocks. Granodiorite, 
granite, quartz diorite. 


*resistivity 5000-20000 2 M. 


Devonian sedimentary and volcanic rocks. 


*resistivity 100-2000 Q M. 


Ordovician and Silurian sedimentary rocks 


AresiStiy) t¥e00-200 900M; 


* The resistivities in regions 1 and 2 were obtained from 


Rankin and Reddy [1969], Peeples [1969] and Cochran and 


Hyndman [1970]. The resistivities for the rocks in other 


regions were obtained from Keller and Frischknecht[1966]. 
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[Caner et al., 1969; Rankin and Reddy, 1969; Peeples, 

1969; Cochran and Hyndman, 1970]. In the mountains of 
southern British Columbia (region 1 on Figure 15) a thin 
(10-20 km) layer of resistive rocks overlies a thick 

(20-40 km) region of rocks with very low resistivities 
(5-15 2M). Conversely over most of southern Alberta (region 
2) a thin (1-2 km) layer of conducting sediments (5-10 9M) 
overlies layers of more resistive (200-400 2M) Cretaceous 
sedimentary rocks [Rankin and Reddy, 1969; Peeples, 1969]. 
There is no doubt that in the nocntain regions induced 
earth currents strongly affect the magnetic fields of Pc 4 
and Pe 5 micropulsations. These effects migh 
far east as Calgary. Induction should also be noticeable in 
the Pc4 micropulsation band at all the stations in southern 
Alberta (see Appendix A7). Highly resistive intrusive rocks 
occupy most of the Canadian Shield (region 3) and, depend- 
Ing on the thickness of these rocks, the induction effects 

at FITCH, SMIT, and RELI should be negligible. Ocean waters 
should produce the most substantial induction at CAMB. There 
is also evidence for a shallow,conducting anomaly to the far 
north on Ellesmere Island [Praus et al., 1971] but the 

exact extent of this anomaly is not yet known. In summary 


we expect induction to be noticeable at all our stations 


except the three on the Canadian Shield. 
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ZoeicenOdaOtsdata analysis 


Many researchers in the field of micropulsations 
estimate wave characteristics such as mean periods, mean 
amplitudes, and polarization hodographs directly from 
amplitude-time recordings. These methods are often cumber- 
some and the results misleading. (See for example Pope 
[1964], Egeland [1965].)In addition, the necessity for 
clean monochromatic events drastically limits the selection 
of data. Finally, the results of such analyses are often 
difficult to compare with the theoretical literature since 
most theories are developed in terms of the spectral com- 
ponents of the wave forms (i.e. the time dependence is of 
the form ae Clearly an alternative presentation of 
the data is needed. 

Correlation and cross spectral matrices offer 
useful, alternate formulations for micropulsation data (see 
Appendix Al). Data pertaining to quasi-monochromatic waves 
are easier to interpret by considering the spectra of the 
waves and the geometric figures traced out by the waves' 
disturbance vector. In the simple case of a monochromatic 
Signal the endpoint of the disturbance vector traces out 
an ellipse. The parameters of the polarization ellipses 


can be computed directly from the cross spectral matrix 


Jig(FyoPy of) (see Appendix Al for an explanation of the 
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symbols) and the relative phases of the vector components 
at different stations can be calculated directly from 
J5 5% %o>f). Considerable care must be exercised in 
using these methods, however, since many hazards are 
involved in estimating the spectral parameters from real 
data (see Appendices A2, A3 and A4). 

In the analysis of the data presented in this 


thesis the following parameters were calculated. 


(A) Polarization Parameters 


(a) R, is the ratio of the apparent plane polarized 
Signal intensity to the total signal intensity 
(see equation (Al17)). The definition of plane 
polarization used here is not the same as that 
used in optics. By plane polarization we mean that 
at one point in space the vector is restricted 
to motion on a given plane. For example, suppose 
that we have sinusoidal pulsations with powers a, 
Demdidecmnlecivern.] 0. and: ZedinectlonsSmrespeGLiVely 
and that a totally random signal, with power d 
in each direction is then superimposed on the 
pulsations. Since the motion of the disturbance 
vector of three-dimensional, sinusoidal oscillations 
describes an ellipse in a given plane, the intensity 


of the plane polarized signal is (at+btc) and 
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ia aed artes) + (ctd i 
Ro is the ratio of the completely coherent signal 
intensity to the plane polarized signal intensity 
(see equation (A20)) R, is equivalent to the 
polarization ratio defined by Born and Wolf [1954] 
and Fowler et al. [1967]. 
The intensity I is the total power of the coherent 
signal (equation (A22)). In the above example 


ip=seat bre: 


The projection of the polarization ellipse (equation(A21)) 


was determined for each of three planes (H-D, H-Z, D-Z) 


by using the rotation matrix R (equation (Al14) and the 


cross spectra Pas (equation (A18)). 


(d) 


In each plane, the ellipticity is the ratio of the 
Minors lOmtieca id) OG axis sOre cle snOlarizatl On enllpse. 
Looking at the plane from a direction in which the 
coordinates form a left handed system (downward on 
the horizontalsH-D plane), the ellipticity is 


considered positive if the polarization is counter- 


mclockwise biguYres 16). Zero ellipticity ~(E=0) 


represents linear polarization while unity represents 


Circular polarization . 
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(e) The polarization angle 6 is the angle, measured 
clockwise, between the major axis of the ellipse 
and the H direction in the H-D plane, the major 
axis and the H direction in the H-Z plane, and 


the major axis and the D direction in the D-Z plane. 


(B) Station to Station Correlations: Coherency and Phase 


(a) The coherency RT gst 12 te! is a measure of the cor- 
relation of component i at station a with 
component j at station ro (see equation (A27)). 
This might be, for example, a correlation of the 
H component at SMIT with the Z component at MCMU. 
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(b) The phase estim 
phase difference between component i at station 
ry and component j at station Po. {Bae Fig is 
positive then component i leads component j. 
To prevent needless repetition in the calculations of the 
phases and coherencies only six component-to-component 
comparisons were made for each pair of stations. These were 
H(r,)-H(r,) H(r,)-D(r5), D(r,)-D(ro), Det, es and 
HY.) 2115) - For the sake of brevity, not all these 
estimates are given in this thesis. In most cases one 
station was used as a base and the coherency and phase 


parameters of all other stations calculated with respect 
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In the actual analysis, two data windows were 
used. One was one hr long and was used in the analysis 
of low frequency Pc 5 micropulsations. The other was 
1/2 hr long and was used in the analysis of Pc 4 micro- 
pulsations. A more complete description of these data 
Windows and the associated spectral windows is given in 
Appendices A2, A3 and A4.3. 

The lengths of the data windows were dictated 
by the scale size of the temporal variations of the 
micropulsation spectral parameters. Normally Pc 4's and 
Pc 5's have stationary spectral characteristics over periods 
Otel/z2e@hour to mours, With perhaps some Fluctuations in 
amplitude over periods of 10-20 min [Satto, 1967]. There 
is generally no consistent and pronounced dispersion in 
the spectra and no consistent short term changes in the 
sense of polarization. The long term stability of these 
parameters is illustrated in Figure 17. Apparently the 
ellipticity and central frequency of this event remained 
stationary for over an hour. -This feature was generally 
the case in almost all the events analysed in this fashion. 
It appears then, that the spectral changes are based on a 
time scale greater than 1-2 hr, perhaps a time scale of 


diurnal length. With these considerations in mind, a window 
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El Gaeta CONCOUYS=: Of Lie mim censa Gy. log19 
(intensity). and ellipticity ine the H=D plane 
Ofscdumicroplilsation event OccurGingson days 163), 


1970 and recorded at LEDU. Intensities are in 


Und tS= orf 2 
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length of 1/2 - 1 hr apparently gives the best compromise 
between temporal resolution and the amounts of data 
analysed. 

Computer programs were written to determine the 
elements of the cross spectral matrix (equations (Al12), 
(A25)) and from these, the polarization, phase and coher- 
ency estimates were computed. Appendix A4 outlines the 
practical considerations in the evaluation of these para- 
meters. To facilitate analysis, the computer programs 
were also designed to pick quasi-monochromatic spectral 
peaks (corresponding to quasi-sinusoidal micropulsation 
events). The quasi-monochromatic restriction specifies 
that the micropulsations have a spectral peak of mean 
crenienty v with a peak width Av (measured between 
half power points) such that a <2 leeninrpractice , 
DeakscewicheeAvyy >) Ol Spewereeres ccrede ey ihesgbesterasolution 
used in separating peaks is given by the relation Av./v > Oral 
where Av. is the separation of the two spectral peaks and 
vy their mean frequency. For example, peaks at 100 and 110 
mHz would be analysed separately wnereas peaks at 100 and 
105 mHz would be averaged to give an estimate at ~ 103 mHz. 
In the selection of the polarization parameters by the compu- 
ter, events with Ry or Ro < 0.75 were rejected (see Appendix 
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estimates, events with all coherencies less than 0.45 were 
rejected. In other words if Kig(tyoro>fy) <a 4 etn en 
the phase estimates in the comparison of station ry with 
station ro at frequency Uy were rejected. This selection 
rule gives an approximate 95 per cent confidence interval 
in the rejection of unwanted random signals. The selected 
polarization parameters were stored on data cards and later 
USedFinsarStatisticallvanal'ysisi-of=thespolarization para- 
Nevers. 

The selection of spectral peaks by the computer did 
not, however, completely alleviate the tedious task of visu- 
ally checking the computer outputs. Often the spectral peaks 
and polarization ratios (Ri, Ro) of an event differed con- 
Siderably from station to station. In extreme cases a 
spectral peak at one station might be nonexistent at 
another station. To find the polarization parameters of 
a given event at all stations, it was then necessary to 
make direct comparisons of the spectra and the amplitude- 
time plots of the micropulsations. 

Only the two horizontal components (H and D) of 
the 1969 data were analysed. Samson et al. [1971] have 
Summarized the results from this data. A more generalized 
analysis was carried out with the 1970 data in that the 


three-dimensional polarization, coherency, and phase para- 
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meters were calculated. We must remember, though, that 
induced earth currents substantially affect the polariz- 
ation characteristics in the H-Z and D-Z planes and that 
horizontal changes in geoelectric properties substantially 
affect the station-to-station phase spectra. These con- 
siderations will be discussed in more detail in the analysis 


of the data. 


2.4 Selection of micropulsation events 


Micropulsation events were selected from periods 


with low to moderate activity (K, < 4) and with evident 


p 
LGRASOGeECMoOmMiLCrODUTSa LION, ACERV ItymaSm ind «| Obean ours 
or more. The Meanook standard magnetograms were used in 
the initial determination of these intervals since the 
sensitivites of the component magnetograms (H = 10.3 y/mn, 
D = 1.6'/mm, Z= 9.4y/mm) allow the detection of most large 
amplitude Pc 4 and Pc 5 micropulsations. Digital plots of 
the data from our stations were then made in the selected 
time intervals. Those micropulsation events which satis- 
fied the following simple criteria were selected for 
analysis. 

(a) There must be at least four stations with no 


data gaps and with timing accurate to within 


~ 5 sec in the interval of interest. 
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(b) There must be apparent micropulsation activity 
at three or more of these stations and this 
activity must last for more than an hour. No 
restrictions were placed on the shape of the 
pulsation waveforms although some events were 
selected specifically because of their sinu- 
soidal appearance. We placed no restrictions 
on the local time of occurrance of the events. 
(c) The micropulsations must have amplitudes of at 
least 2y (peak to peak on one component) at 
three or more stations. Lower amplitudes lead 
to sizeable quantizing errors in the spectral 
estimates (see Appendix A5). Due to the natural 
falloff of power in the geomagnetic spectrum this 
restriction limited the highest selected frequencies 
COg Varsunmi 2s 
Table 3 lists the events selected for analysis, and Table 
4 gives a summary of the KO indices during the selected days. 
Indices corresponding to the time intervals analysed are 


underlined. 
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Indices on the Days Selected for Analysis 
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3. OBSERVATIONS 
Sete Ur OUuuUC.c On 


The data presented here occupy two distinct 
SeCulolsmmm hemi icStasec tion .cOnballicmanrs cacis ticad | 
analysis of the three-dimensional polarization parameters 
of the selected micropulsation events (see Table 3). The 
second section contains a complete, detailed discussion 
of a number of individual micropulsation events. The 
statistical results are presented first in order to develop 
a good perspective for the analysis of the complicated 
characteristics of individual events. 

Tne statistical anaiysis summarizes tne average 
characteristics of the intensities, ellipticities and 
polarization angles.(see Section 2.3) of quasi-monochromatic 
Pe 4's and Pe 5's. The variation of these parameters with 
time, latitude, and frequency is considered. Since the number 
of parameters is very large we will not attempt to discuss 
all the data but will present only those results which are 
representative and illuminating. The selection of polar- 
ization parameters is discussed IN edetalM ioe ClLons css 
and A4 and will not be repeated here. Unless otherwise 


indicated all the statistical results pertain to the 1970 


data. 
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Figure 18 summarizes the total number of data 
hours in each hour of Universal Time. Histograms (1) and 
(2) pertain to data which were subjected to polarization 
analysis using the 1/2 hour and 1 hour data window res- 
pectively. Histograms (3) and (4) pertain to data which 
were subjected to station-to-station phase comparisons 
using the 1/2 hour and 1 hour window respectively. Most 
of the data are found during the daytime hours reflecting 
the local-time-dependence of the occurrance frequency of 
the micropulsations. Normally data from the morning events 
are the clearest and the easiest to analyse, showing very 
predictable ellipticities and polarization angles at all 
the stations. Because of the relatively small amounts of 
data from the nighttime, inferences from data during this 
time should be considered very cautiously. There is also 
the problem that some of the nighttime events are composed 
of short-lived impulsive trains much like Pi's. These 
events may have characteristics that are fundamentally 
Gitirerent surom -Cu4mom PC Se characteristics. 

Data from MENK are included in the analysis even 
though the gain of the D component at this station changed 
sporadically. To correct this varying gain we compared 


our D component data with standard magnetograms from Meanook. 
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After calculating an average value for the gain of the 

D component we corrected the MENK data by multiplying 

by this average gain. On an event to event basis this 
leads to rather variable characteristics. The statis- 
tical averaging of the events should, however, give mean- 
ingful results. In fact we find that the sense of polar- 
UZetiOnmd etl S ss taclOlme 1s eCONSISUtCIC wit suneeresultSeat 
the other stations although there are some variations in 
intensities, ellipticities, and polarization angles 
calculated at this station. 

A complete discussion of the polarization and 
station-to-station phase characteristics of five represent- 
ative events is given after the statistical analysis. The 
frequencies of the micropulsations in these events range 
from 1.3 mHz to 20.3 mHz and thus cover the whole spectrum 
of interest here. The events all occurred during the day- 
time, the earliest between 1400-1700 UT and the latest 
between 2300-2400 UT. The event on day 186 is the most 
intriguing, having very sinusoidal oscillations with peak 
to peak amplitudes of nearly 180 y in the H component at 
SMIT. These large amplitudes permit the use of this event 
in a worldwide analysis of the spectral characteristics. 
The longitudinal dimension inherent in a worldwide analysis 


is especially helpful in evaluating the statistical, diurnal 
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characteristics obtained from the data at our line of 
stations. We accordingly digitized magnetograms from 
13 other magnetic observatories in the northern hemisphere 
and used this data to determine the polarization character- 
istics of the event at those observatories. Since standard 
magnetograms normally have poor temporal resolution, these 
polarization characteristics must be considered with care 
In the presentation of the data that follows 
corrected geomagnetic coordinates are used throughout. All 
latitudes in the text and all latitudes in the figures are 
corrected geomagnetic latitudes. Hakura (1965) gives a 
complete description of the computation of these coordinates 
on a worldwide basis. For simplicity we have plotted all 
temporal variations in Universal Time. Unless otherwise 
indicated, Universal Time is used throughout the text. An 
approximate conversion to local geomagnetic time (LGT) at 
our stations can be made by substracting 8 1/2 hours from 


the given Universal Times. 


SecemAaS tatisticalvanalysisvofethetpolarizati onspanameters 


In the analysis of the data it soon became apparent 
that almost all Pc 4's and Pc 5's have frequencies which are 
latitude-independent. No gradual latitudinal changes in 


frequency are evident, though two or more micropulsation 
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trains of different frequencies often occur at the same 
time. This feature is evident from a visual inspection 

of the magnetograms, from power spectra calculated for 

the stations, and from the spectral peaks selected through 
use of the computer programs (see Section 2.3). This 
evidence supports the observations of Fllis [1960] and 
Obertz and Raspopov [1968] who also found no latitudinal 
changes in frequency. Exceptions to the above rule do 
occur, however, and an example of an event with latitude- 


dependent frequencies is given in Section 3.3. 


SCS I Ntens loves mie miSOmIncenSicyecOntOUns «ini lgure 
19 summarize the daily variations in the relative inten- 
ShCLese0 feel Seniiize MLCrOopUlus dil ONSmOCCU Tat NGOs 11) OUTELO sa 
data. We have used only those events with distinct 
spectral peaks (selected by the computer) at three or 
more stations. The intensities at the other stations were 
determined by visually inspecting the computer outputs 

of the spectral parameters (see Section 2.3). The 
intensities at each station were then summed in hourly 
intervals and the average intensities, Is, in these 
intervals calculated. If L is the largest of the 
intensities at all the stations in a given interval, then 


the relative intensity at station i is defined by I,/I.. 
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The intensity is clearly localized in latitude and moves 
northward during the day. The maximum intensity is at 
65°-68°N during the early morning, moving to ~u 75°N 
between 0830 and 1530 LGT (1700-2400 UT). These diurnal 
variations closely resemble those of the auroral oval. 
The resolution in the northern regions is poor since CAMB 
and SMIT are 9° apart and consequently the position of 
the maximum is uncertain in the interval 1700-2400. 
Equivalent diurnal trends also appear in the intensity 
maxima Of 3-8 mHz Pc 5's but the maxima are south of 


thoseeinerigure 19>" Pc 4 micropulsations, on the other 


hand hava wv 
seKnremy ee on on ry 


ather complex iso-intensity contours. 

Latitude plots of the relative intensities of 
micropulsations in given spectral bands indicate that 
low frequency Pc 5's have most of their energy concentrated 
north of 70°N, whereas Pc 4's have their largest intensities 
and most energy south of 60-62°N. Figure 20 presents 
latitude plots of relative intensities in 7 spectral bands 
ranging from 0-2 mHz to 15-20 mHz. To obtain these plots 
we calculated the relative intensities for each individual 
event occurring in 1970, sorted the events into the 7 


spectral bands, and then averaged the relative intensities 


at each station. We have not included MENK data in this 
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pulsations in seven spectral bands. The shaded areas 
are the estimated positions of the intensity peaks. 
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plot because of the anomalous D component at this station. 
The shaded areas in Figure 20 outline the regions where 
we estimate the peak intensities occur. In the 0-2 mHz 
spectral band the shaded region is a southern limit and 
in the three Pc 4 spectral bands the shaded regions are 
northern limits. Clearly the latitude of the intensity 
peaks depends on the frequencies of the micropulsations, 
the peaks for the highest frequencies occurring at the 
southernmost latitudes. One difficulty in interpreting 
the data in these plots arises from the fact that the 
averages are taken from events occurring at any time 
during the day. Diurnal effects such as those noted in 
Figure 19 definitely broaden the latitude ranges of the 
intensity peaks and make the peaks more difficult to 
resolve. 

An alternative method of depicting the intensity 
data is to construct scatter plots showing the estimated 
latitudes of the intensity peaks of micropulsations of all 
frequencies. Figures 21 and 22 show the results of such an 
anaiysis. In making these scatter plots the day was 
divided into two sectors, based on the diurnal variations 
noted in Figure 19.The first sector extends from 0200 to 
1600 and the second from 1600 to 0200. We determined 


the latitudes of the intensity peaks by noting the station 
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Fig. 21. Estimated latitudes of the intensity peaks 
of micropulsations occurring between 0200 and 1600 UT. 
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Fig. 22. Estimated latitudes of the intensity peaks 
of micropulsations occurring between 1600 and 0200 UT. 
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with an obvious maximum in intensity for a given micro- 
pulsation event. If the intensity showed peaks near two 
different stations then the latitudes of both these peaks 
were plotted, even though one peak might be larger than 
the other . If there were more than two peaks or if there 
were no clear peaks, we rejected the event. Since the 
two northern stations are separated by 9° the latitudes 
of peaks occurring between these stations were determined 
by extrapolation, the peak being proportionally closer 

to the station with the greatest intensity. In some 
instances the intensity increases monotonically across 
the array of stations. YJhese events, which clearly have 
maxima beyond the latitude of the array,have their peaks 
plotted at 58°N or 77°N depending on the direction of 
increasing intensity. 

The most obvious feature in each plot is a strong 
latitude-frequency correlation in the grouping of the 
estimates below 72°. In Figure 21 (0200-1600) the peaks 
are centred near 10 mHz at 60°N, decreasing continuously 
to about 2 mHz at 71°N. In Figure 22 (1600-0200) the 
peaks are centred near 15 mHz at 60°N, decreasing contin- 
UOUST¥ato 2 mz at 70-75°N. A cluster of intensity peaks 
which covers all frequencies and is localized in northern 
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region corresponds to the approximate position of the auroral 
OValSduUrrigethus timnewminteyval CUseceSechiant4) i). ‘alheshigh 
Tatitudercluster does not appear in Figure 21 but this might 

be expected since the nightside segment of the auroral oval 

1S USUCMLY Eat bbatitudes south of 70°N. The absence of any high 
latitude cluster in Figure 21 also indicates that nightside 
micropulsations might have different energy sources than those 


Causing dayside micropulsations. 


Ue eee TOLLE Les Figures 23 to 29 summarize the diurnal 
Grinders tressUh PCeseand hces Gliipu.creies.. | MOS t Of tire 
data are from 1970 although some of the H-D plane data are 
from 1969. The estimates used in this analysis and in the 
analysis of the polarization angles were all selected by com- 
DUCES TeCctmtoO tne UCOnStraints 1isved) inesectLon 2.3. einese 
plots are based on averages over four-hour intervals and are 
Calculated tor ali 3 planes (H-D, H-Z, D-Z)2 [he circles in 
the plots are at the centres of the four-hour intervals. 

The ellipticities at MENK in the H-D and H-Z planes 
appear somewhat more linear than those at LEDU and CALG. This 
feature is probably due to the malfunctioning of the D compo- 
nent at MENK. The diurnal trends in the ellipticities at MENK 
do, however, agree with those at CALG. and LEDU. 

In the H-D plane, the low frequency Pc 5 micropul- 
sations (0-6 mHz) show the simplest and most consistent po- 


larization patterns. At all stations from SMIT southward 
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Fig. 23. 


A summary of the ellipticities in the 


three planes of micropulsations recorded at CALG. 


The data points are based on averages in 4-hr 


intervals. 
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the sense of, polarization is, CC. (counterclockwise) over 
most of the day with a period of CW (clockwise) polar- 
ization occurring between 2000-0400 (1130-1930 LGT). 
Tne period of CW polarization is most, evident at FICH, 
WNeremtiere lI pti Cityedecy.eaSeSe t0s- 0-5. wand Ss. 10 lar 
ization pattern corresponds remarkably well to those 
computed by Kato and Utsumi [1964] and Rankin and Kurtz 
[1970]. The micropulsations at CAMB are the exception 
to this diurnal polarization rule, having CW polarization 
from 1200-2000 and CC polarization outside this time 
Hitervdleersioese diurnal. polarization characters s ics 
will be discussed in greater detail later. 

Pc 4 micropulsations, on the other hand, appear 
to have rather complex polarizations in the H-D plane. It 
should be noted that the data used to construct the Pc 4 
plots are rather limited, especially during the nighttime 
hours (0200-1400). Morning Pc 4 micropulsations appear to 
have some predictable polarization characteristics and these 
Nishigdisoebemdiscussed in detail Vater. 

The low frequency Pc 5's show relatively simple 
polarization changes in the H-Z plane. Micropulsations 
at most of the stations have CW polarization over the whole 


day, with short periods of CC polarization occurring in 
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the interval 0800-1200 ate CAMERON cond ra GH a de 
patterns at MCMU, MENK, and LEDU are the simplest having 
almost constant CW polarization throughout the day. With 
tihemexCepeloneot a short periodeotmcGrepolarization) 1m tne 
O-8 miz band from 0400-1200, (Pe 5 “sat CALG al'so have 

CW polarization throughout the day. The polarization 
patterns at CAMB are somewhat more complex than those at 
the other stations but distinct similarities between CAMB, 
SMIT, and FITCH do exist, especially in the 0-4 mHz band. 

In the H-Z plane changes in the ellipticities 
of Pc 4's have more complicated characteristics than do 
thosemoterc.o Ss.) ihe Pc 4 *patternseares also coniused by 
sketchy data in this spectral band, especially at the 
northern stations. Generally speaking, polarizations 
at the southern stations are very linear, with little 
apparent structure. This characteristic may be a mani- 
festation of earth induction effects in the southern 
prairie regions. 

Polarizationsin the D-Z plane show the least 
predictable characteristics. At CALG and LEDU polarizations 
in the Pc 5 band change slowly over the day. At these two 
stations the greater part of the day has CC polarization, 
with the exception of a short period of CW polarization 


at LEDU in the interval 0400-1200. Polarizations are 
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most elliptical around 2000-2400 at which time the 
Ci piolctlLy ese tO. 3 ateCALGeandeleDUm The: perivodlot 
CW polarization becomes more pronounced at MENK and 
extends, from 0800-2000. At SMIT sands Fach, Pe 5 polan- 
ization patterns are very complex, changing rapidly 
with time and frequency. The pattern in the 0-4 mHz band 
at SMIT resembles that in the H-D plane. More precisely 
polarizations are CC throughout most of the day with 
a short section of CW polarization from 2000-0400. Pc 5 
micropulsations at CAMB show predominantly CW polarization 
with a short section of CC polarization occurring between 
1600 and 2000. Once again micropulsations in the Pc 4 
band show predominantly linear polarizations at the 
southern stations. The linear polarizations are almost 
always due to the small Z components at these Se eS, 
Figure 30 illustrates the temporal polarization 
features of the very low frequency Pc 5 micropulsations 
(0-4 mHz). The samples are based on averages over two- 
hour intervals and the circles in the plot are at the 
centres of these intervals. This figure shows that these 
micropulsations obey rather simple polarization rules in 
all three planes. In the H-D plane the ellipticities 


follow an inverted U-shaped diurnal variation (based on 
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a universal day) at all stations except CAMB. Polarizations 
are most elliptical in the period 1000-1800 having ellip- 
GiGi comdSmmlargemas e+ oma tao ieed lant CHem Av dil iss ce. b1-0N's 
from SMIT southward micropulsations are polarized in a CC 
sense during the evening and morning hours with a 7 to 8- 
hour period of linear and CW polarization occurring between 
2000 and 0400 (1130-1930 LGT). This corresponds to the 
period of CW polarization noted by Kato and Utsumi [1964] 
fOGsh Ceo Smatenlone lati cuceas tabi Ons. ss UUrINGuthe Gay tiie, 
micropulsations at CAMB have a reversed sense of polariz- 
ation as compared to the other stations. This implies that 
during the daytime there exists a given latitude across 
which the sense of polarization reverses. We shall refer 
to this latitude as the demarcation line. The existence 
of a demarcation line in this frequency band suggests 
that we plot latitude profiles of ellipticities in other 
frequency bands to determine whether this reversal exists 
for all frequencies. This problem will be dealt with in 
detail later. 

Polarizations in the H-Z plane are predominantly 
CW throughout the day. Pulsations at stations from MCMU 
southward always remain CW whereas pulsations at FITCH, SMIT 


and CAMB show periods of pronounced CC polarizations at 
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0600-1000, 0600-0800 and 1000-1400 respectively. There 
is no evidence for a consistent latitudinal reversal in 
the sense of polarization as there is on the H-D plane. 

Ellipticities in the D-Z plane are somewhat 
more complex than those in the other two planes but some 
predictable patterns are evident. At CALG, LEDU, MENK 
and MCMU the ellipticities show a U-shaped diurnal 
variation, being mostly CC from 2200-0200 and CW from 
0600-1400. At FITCH, SMIT and CAMB the polarizations are 
more confused. Micropulsatons at FITCH are linearly 
polarized for most of the day, whereas at SMIT they show 
a diurnal variation very similar to that in the H-D plane. 
In fact, Pc 5's at both SMIT and CAMB have diurnal vari- 
ations which are the reverse of those at the southern 
stations. More precisely they have CC polarization in 
the morning hours (0600-2000) with a tendency toward 
CW polarization in the afternoon and evening. 

Figure’ 3 Shows, plots OF ‘ellipticity versus 
latitude in the H-D plane for the morning hours 1200-1600 
and 1600-2000. Data from MENK are not included since the 
CIM DLACIUTeS -duminatsSitacvon aresinaccumate. @ac S micro- 
pulsations show an evident reversal in the sense of 
polarization occurring somewhere between 70°N and 75°N. 


Unfortunately the station spacing in this region is very 
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poor (CAMB and SMIT are 9° apart) and consequently we 
cannot determine from these plots if there is any latitude- 
Frequency dependence of the demarcation line. Obertz and 
Raspopov [1968] noted that the demarcation line apparently 
moves southward as the frequency of the micropulsations 
increases. Samson et al. [1971] have also shown that for 
selected, individual events there appears to be a latitude- 
frequency dependence of the demarcation line. In the Pc 4 
band the latitudinal characteristics of the ellipticities 
appear more complex. There are, however, indications 
that these micropulsations have a demarcation line between 
60° and 65°. This feature is most evident in the interval 
1200-1600. 

We can obtain additional, useful information 
DVS LObCINdaGInTULUChLY aS eas GUNCLLONGOh tT requency «for 
the seven stations. Figure 32 indicates that a frequency- 
dependent change in the sense of polarization in the H-D 
Dlanew1s evident cat SMIq” FICHS and MCMU.s) At SMIT micro- 
pulsations with frequencies less than 7 mHz have CC 
polarization while those with frequencies greater than / 
mHz have CN polarization. At FICH and MCMU the change occurs 
at ~ 9 mHz. CAMB does not show any distinct changes in 


polarization and micropulsations at most frequencies 
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Fig. 32. Ellipticities in the H-D plane plotted 
against frequency at 7 stations. All the data are 
from the interval 1200-1600 UT. The solid lines 
have been fitted to the data by taking averages 


over 2 mHz bandwidths at intervals of 1 mHz. 
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are nearly linearly polarized. There is, however, some 
indication of a change at 2-3 mHz. The average ellipt- 
icities at the three southern stations (MENK, LEDU, and 
CALG) do not show any change in sign. If a change does 
occur at these stations it must be only at frequencies 
greater than 15 mHz. In fact, Rankin and Kurtz [1970] 
found that at Leduc a change in the sense of polarization 
occurs near 20-40 mHz. Some caution should be used in 
interpreting the Calgary data, however, since there are 
a considerable number of events above 9 mHz which have 
GHevolarization. gel lipticity plots aromedataerecorded 
in the late mornina and afternoon indicate that the change 
in sense is seen at higher frequencies at this time (see 
Figure 33). This trend might be expected from the diurnal 
variation of the intensity maxima shown in Figure 19. 
LheapLlots.inkigures32-alsoweus-.toses timate 
the latitude of the polarization demarcation line as a 
function of frequency. The frequency at which the sense 
of polarization changes (e.g. 7 mHz at SMIT) corresponds 
to the frequency of micropulsations that have a demarcation 
line at that station. These estimates eliminate the poor 
latitude resolution inherent in ellipticity-latitude 


profiles computed from data from widely separated stations. 


The weakness of this method is that it demands large 
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amounts of micropulsation activity at one time or station- 
ary geomagnetic conditions over long intervals. Neither 
condition prevails very often and consequently the data 

are scattered and frequency resolution is lost. Figure 

34 shows the estimated frequencies of the polarization 
reversals plotted against latitude. The solid line gives 
the approximate trend of the intensity peaks plotted in 
Figures 21 and 22. Evidently there is a correlation between 


the ellipticity reversals and the intensity peaks. 


Sweet Olarization ang )es A simple plot of average 


polarization angles is not possible because of the cyclical 
nature of these parameters. The distributions of the angles 
are conveniently presented in the form of circular bar 
histograms. The lengths of the bars give the expectation 
of finding the polarization angle in the given intervals. 
For the data presented here we have used the six intervals 
SHOWN Lin Eb POULe @oOeemen Dar oOte Lendcimlavuernducates® that 
all measured angles were in the given interval. For clarity 
in the following histograms bars with lengths of less than 
Onalmare notaploutved™ 

Not all the data will be presented here since 
this would require the construction of over 1200 histograms, 


based on intervals of 4 hr x 2 mHz. In general the polar- 
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Fig. 34. Estimated frequencies of the polarization 
reversals at the stations. The arrows at CALG (58.7) 

and MENK (62.59) indicate that the reversals occurred 

at frequencies higher tnan those marked by the horizontal 
bars. The frequency range of the reversals at LEDU (61.29) 
was determined from the data given by Rankin and Kurtz 
[1970]. The solid line shows the approximate trend of 

the intensity peaks plotted in Figure 22. 
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Fig. 35. The ranges of the polarization angles in 
the circular histograms. In tnis example we have 
assumed that all the polarization angles are in the 
0-309 range. 
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ization angles of the low frequency Pc 5's are the most 
predictable; the angles in the Pc 4 band show considerable 
SCocceoremeiethesh- le pianes thes pollamizavnon angles 01 PC. > 4s 
indicate that the magnetic fields at CALG are strongly 
affected by induced earth currents, the most probable 
cause being a deep conducting anomaly in the Cordillera 
region. The following histograms give summaries of the 
H-D, H-Z, and D-Z polarization angles of 0-4 mHz pul- 
Sations and a summary of the H-D polarization angles of 
micropulsations in the 10-15 mHz band. We have omitted 
the Pc 4 data in the H-Z and D-Z planes because of the 
small amplitudes of the Z component in this band and 
because of induction effects on the Z components. 

In the 0-4 band, the polarization angles in 
the H-D plane (Figure 36) show simple diurnal variations 
at all stations except CALG. The H direction predominates 
over most of the day with tendencies for pvolarization in 
the D direction at CAMB, SMIT, and FITCH between the hours 
of 0800 to 2000. This interval marks the time during 
which pulsations at these stations have the most elliptical 
polarization in the H-D plane (Figure 30). The tendency 


forscobarization inethe D dinection™is most pronounced in 


the interval 0800-1200 at CAMB changing to 1200-2000 at 
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SMIT and FTCH. The major axes at MCMU and LEDU are 
clustered around the H direction throughout the day and 
show only slight diurnal variations in direction. From 
0400 to 1200 the angles favor the 0°-30° interval and 
for the rest of the day they favor the 150°-180° interval. 

The polarization angles at CALG differ sub- 
Stancuial by 1 romethose watethe other.stavions. lhe wmajor 
axes are grouped in the 1a Sse region and there is 
little variation over the entire day. In fact, the major 
axes are parallel to the strike dividing the plains from 
the mountains. The direction of the strike is v 30° west 
of geographic north or 53° west of magnetic north (cor- 
responding to a polarization angle of NEPAL. This type 
of polarization would be expected if CALG were on the 
resistive side of a vertical geoelectric discontinuity 
(see Appendix A7, Figure A21). This supposition is 
supported by magnetotelluric and geomagnetic depth sounding 
datauinesouthern: dA libentay(see Section. 2.c). 

The polarization angles in the H-Z plane (Figure 
37) have very simple diurnal characteristics. At all 
stations the preference is for the H direction although 
there is some change of angle with latitude. At FITCH, 


MCMU, LEDU and CALG most angles lie in the 0°-30° interval 
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whereas at SMIT the angles are in the 150°-180° interval. 
These histograms and the ellipticity plots show that the 
pulsation magnetic fields are definitely not polarized 
transverse to the main geomagnetic field. For example, 
transversely polarized waves at LEDU would have polariz- 
ation angles grouped around 167° whereas the actual angles 
are grouped around 15°. Pc 5's at CAMB show more confused 
characte risi1 css thane those atthe other Stations. “Even so 
the major axes of the polarization ellipses generally point 
ineune He direction. 

The characteristics of the polarization angles 
in the D-Z plane (Figure 38) show pronounced diurnal 
variations at every station except CALG. Conditions in 
the interval from 1200 to 1600 are the clearest with the 
magormaxes at most Stations favoring the D direction. 
Micropulsations occurring during this time interval must 
have large D components since they are also polarized in 
the D direction in the H-D plane. At other times of the 
day, the polarization in the D-Z plane is more confused, 
especially at FITCH and MCMU. During these times the D 
and Z components are probably of comparable magnitude. 

The intervals of confused polarization are most pronounced 


between 0000-0800 at FTCH, MCMU, and LEDU. The polar- 
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Fig. 36. Polarization angles in the H-D plane 
of micropulsations in the 0-4 mHz spectral band. 
The sample intervals are 4 hr long. The horizontal 


Ddteatetie top of the diagramelssieunit.11ong. 


Fig. 37. Polarization angles in the H-Z plane of 


micropulsations in the 0-4 mHz spectral band. 
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lize meal Guiesm 0:1 ah Clo.. Seno te LeCumcls Ona Vy Once 
direction in this time interval. On the other hand, 
thepanglesmof Pc 5's at CAMB favor the Z direction 
in the interval 2000-2400. 

Figure 39 summarizes the diurnal polarization 
characteristics in the H-D plane of low frequency (0-4 
mHz) Pc 5 micropulsations. Data from CALG are not 
included in this plot because of their anomalous behavior. 
The dotted line on the plot gives an estimate of the 
position of the intensity maximum as a function of time 
(see Figure 19). The polarization patterns show some 
strikingly simple features, espectally when they are 
correlated with the diurnal variation of the intensity 
maxima. The major axes of all the ellipses tend to 
line up perpendicular to this line. This is most evident 
for those stations closest to the intensity maxima (CAMB, 
SMIT, and FTCH) although the trend does show at MCMU 
and LEDU as well. Generally speaking, polarizations become 
most elliptical when a given station is near the line of 
maximum intensity, although there is also a tendency for 
more elliptical polarizations at the southern stations in 


the morning hours (0800-2000). In addition, if the line 
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0-4 48 8-12 [216 16-20 20-24 
TIME (UT) 


Fig. 39. A summary of the polarization parameters 
in the H-D plane of micropulsations in the 0-4 mHz 
Snveéctralsband. LELIIpses With solids lanesedepict CC 
polarization and ellipses with dotted lines depict 
CW polarization. The dotted curve shows the position 
of the intensity maximum. The sizes of the ellipses 
are not scaled to the average intensities and all the 
major axes have the same length. The centres of the 
ellipses are at the latitudes of the stations from 
which the data were obtained (CAMB, SMIT, FITCH, MCMU, 
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of maximum intensity is midway between CAMB and SMIT, 
both SMIT and FTCH have a sense of polarization opposite 
to that at CAMB. In the early morning hours, when the 
aGtivicy centre 1s) nears SMim (0400-1200 )mal mistations 
including CAMB have the same sense of polarization. Some 
caution should be exercised in correlating the polarization 
angles of Figure 39 with the line of maximum intensity 
S¥ynce the polarization angles inithe figune=are spatial 
Orientations whereas the intensity line is a plot of 
temporal variation. 

We conclude this study of polarization angles 
Wath a Snort discusstoir of ‘the polarizations “of Pc 4"s in 
the H-D plane (see Figure 40). The characteristics of 
these micropulsations are rather complex and our data 
sketchy so that we will mention only a few of the more 
evident trends. The most noticeable feature in Figure 40 
TSacne CONnSistent polarizations ins thesDedirection, at CALG: 
This feature may once again be a manifestation of induced 
earth currents, although the angles do not align with the 
strike of the Rocky Mountains. LEDU mile Shows a predom- 
inance of angles in the D direction so it may be that these 
DOlamizacOns aremtruly dereimecclOneOi ms ulemoourCenCOnimhg — 
urations. All the other stations have very complex 


polarization characteristics, which differ from station to 
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Fig. 40. Polarization angles in the H-D plane of 
micropulsations in the 10-15 mHz spectral band (Pc 4). 
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Station and from time to time. The polarization angles 

do, however, favor the H direction during the interval 
0000-1200 and the D direction during the interval 1200-2400. 
Pc 4's at FITCH show this feature most clearly whereas at 
MCMU the angles tend to align with the H direction through- 


out the day. 


3.3 Some individual micropulsation events 


This section contains the analysis of five events 
which were individually selected for their clear sinusoidal 
micropulsations. The following pages outline the general 
spectral characteristics of these micropulsations, includ- 
ing the power spectra, station-to-station phase spectra, 
and polarization parameters. 

To reduce the amount of data presented in this 
section the only phases we have plotted are those given 
by F.5(%) sfosf) (see Section 2.3). ~CAMB™15 Used as the 
reference station and all phases are calculated with respect 
to this station. The estimates Faq: Foo and F 33 compare 
the phases of the H, D and Z components respectively at 
Aibestations. “For example, 17 the H components at ETCH 
and SMIT lead the H component at CAMB by 20° and 40° 
respectively, then F,, at FTCH is 220°} "and Po cpcken sli 
fs)-40°. In plotting the phases, we have subtracted the 
corresponding median value from each estimate of F,,, Foo. 


and F334 in order to centre the data on the grids. Thus 
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themes timates oT Fe. at CAMB often do not have the value 
ae. 

It must also be remembered that the phases of 
the magnetic fields can be severly affected by horizontal 
changes in the geoelectric properties of the earth (see 
Appendix A7). Large latitudinal phase shifts in the data 
may not be due to the source configuration of the micro- 
pulsations but may perhaps be due to a vertical contact 
separating two regions with very different conductivities. 
However, we did not note any pronounced and consistent 
phase changes in our data at any given latitude. The 
latitudes of the most apparent phase changes varied from 
event to event, and probably reflected the source fields 
Of the micropulsations - 

The events selected for analysis here are from 
five separate days in 1970 (days 163, 167, 186, 195 and 
262). All the events occurred during the daytime, the 
earliest event occurring at 1400 (0630 LGT) and the latest 
epee oOG I4s05 GI). 

As already mentioned, the frequency of a given 
MPorOUuSalion: eran ts Usual Ny independent Of mic ait e.e seu 
BiSSeXCOOCLON Oftic GVient, OCCURM ING Onmucyecocs tills 1S tne 
case for all the events discussed here. Analysis of the event 
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component apparently decreases with increasing latitude. 
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The event on day 186 merits special attention 
because of its sinusoidal appearance and large amplitudes. 
The magnetograms from 13 standard observatories were 
digitized and these data are included in the discussion 
Oleinismevent. eo ltets intenes tingstomcomparestherspatial 
features of this event with the average polarization 
Characterastics: of..the 19690 and. 1970edata (Figure, 39). 

The events are discussed in chronological order 
and the magnetograms are reproduced in the following 


pages (Figures 41 to 45). 


(a) Day 163, 1600-1700 UT 


This=interval shows distinct 12 mHz,Pc 4 micro- 
DUlsationssanc Lars tinct 1.0 mHz,Pc Semicrcopulsations: 
The Pc 5 activity is most evident at CAMB where it attains 
peak to peak amplitudes of ~v 30y in the D and Z components. 
Note that the D component at MENK may be too small because 
of the malfunctioning of the magnetometer. The Pc 4 activity 
has a distinct latitude-dependence in the powers in the 
three components (Figure 46A), especially in the D component. 
The D components at LEDU and CALG are far larger than those 
ateanveotier Of tne stations.. Converselymtneseceovevent 
has its greatest power at CAMB (Figure 47A), with very 
little power at the southern stations. The H-D ellipticities 
of both micropulsation trains show marked latitude reversals 


(Figures 46C, 47C). The change in the sense of polarization 
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Fig. 41. Magnetograms of the event occurring on 
day 163. The records have been detrended with a 
F=205MHZ digi tel filcer.(see sectionen4. 3)... eile 


time scale is marked in hours of Universal Time. 


Fig. 42. Magnetograms of the event occurring on 
day 167. The records have been detrended with a 


V=20emHz filter. 


Fig. 43. Magnetograms of the event occurring on 
day 186. The records have been detrended with a 


J=207 mz Er iter. 


Pig. 44. MagnetogruansmOiethemevents OCcCUrEINg on 
day 195. The records have been detrended with a 


2eo-100=mMHzZe fT isicer- 


Fig. 45. Magnetograms of the event occurring on 
day e202. mine records have Deen de brenced Ww ithiea 


1-20 mHz filter. 
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Fig. 46. Latitude-dependence of the spectral parameters 
Ciatiewlzemizs mi cropwisati onsm( day slo3ne 1000-1630). 8 ie 
phases of the Z components are not plotted because of 
the low amplitudes of this component at the southern 
stations. 


(A) Power plotted as a function of latitude 
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(B) Relative phases of each component 
@ = Fay (H component phases). 


i => (D component phases). 
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A = F3, (Z component phases). 


The estimates Fiy> Foo» and F33 compare the 
phases of the H, D, and Z components respectively. In 
plotting these estimates we have subtracted the cor- 
responding median value from each estimate of Poe in 


order to centre the data on the grid. 
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(D) Polarization angles. 
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of the Pc 4 event occurs at 63-65°N whereas the change 

in the sense of the Pc 5 event occurs between SMIT and 
CAMB (68-77°N). Inspection of Figure 46B indicates 

that the polarization reversal of the Pc 4 event is caused 
by a rapid change in the phase of the H component, occur- 
ring between MCMU and LEDU. A more abrupt phase change 
occurs between SMIT and CAMB but both the H and D phases 
change by the same amount. Consequently the sense of 
polarization is conserved. On the other hand, the phases 
of the H and D components of the Pc 5 event remain constant 
at the southern stations (Figure 478). Between SMIT and 
CAMB there are sudden changes in the phases of both 
components but the changes are opposite in sign and con- 
sequently the sense of polarization at CAMB differs from 


thacede tae, ovner Stations. 


(b) Day 167, 2300-2400 UT 

Thise1Sea very.simple event, consistingeof only 
4 or 5 sinusoidal pulsations with a mean frequency of 
4.5 mHz. The pulsations are most evident in the H and Z 
components. A plot of power versus latitude (Figure 48A) 
shows that the centre of activity is sharply defined with 
its maximum at FTCH in both the H and Z components. At 


all the stations except CAMB and CALG the power in the D 
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component is much less than that in either the H or Z 
components. This observation is consistent with the 
statistical results for the polarization parameters in 
inatawesexpect mMicropulsations occurring. in this time 
interval to be nearly linearly polarized in the H 
direction on the H-D plane. The pulsations at the south- 
ern stations have CC polarization (Figure 48C) even 
though we might expect CW polarization from the statistical 
GOSULCS eel Nereality, ctoougms the: aftemnoonasectoum 1s 
complicated, showing a large number of events with either 
sense Of polarization. When many events are averaged, 
however, a slight tendency for CW polarization is evident 
(see Figure 30). This event has a change in the sense 
OfepOhahizacion,tnythe H-D planes soccurringsnear GICH 

or SMIT. The H-Z and D-Z planes also show changes in the 
sense of polarization, with pulsations at CAMB having 
DOLdeizachons eOppoOsitestosthose atathessouchernestataons:. 
Figure 48B shows that both the H and D components have 
substantial changes in phase between SMIT and FITCH but 


in this case it is the D component which has the largest 


change. 
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(o)eDay 186 1400-1700 SUT 


This event has an unusual series of sinusoidal 
oscillations with peak to peak amplitudes of 180y in the 
H component at SMIT. The large amplitudes and low frequency 
of this event (= 2 mHz) allow the use of standard magneto- 
grams in the evaluation of the morphological features of 
the pulsations. Consequently we decided to analyse this 
event on a worldwide basis and obtained standard magneto- 
grams from 16 stations around the world. At 13 of these 
stations, the temporal resolution was accurate enough, 
and the pulsations large enough to allow fairly precise 
spectral analysis of the data. These stations are listed 
in Table 5 and their positions mapped in Figure 49. We 
subsequently digitized the magnetograms at 1 min intervals 
and used the digitized data in an analysis of the polar- 
ization parameters of the micropulsations. Magnetograms 
plotted from the digitized data are shown in Figure 50. 
No attempt has been made to ascertain the true accuracy 
of the relative timing of the components at these obser- 
vatories so that the data must be evaluated with caution. 
It is interesting, nevertheless, to compare the Jongi- 
tudinal characteristics of this event with the diurnal 


averages of the polarization parameters (Figure 39). 
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TABLE 5. Magnetic Observatories used in the Analysis 


of Data for Day 186 


Corrected 
Observatory Code Name Geomagnetic Coordinates 


Latitude Longitude 


College 
Siuka 
Victoria 
Newport 
Resolute 
Baker Lake 
Council 


Great Whale 
River 


Ottawa 
St. John 
Leirvogur 
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Fig. 49. The positions of the observatories used in the 
analysis of the micropulsations occurring on day 186. The 
circular grid lines show the corrected geomagnetic lati- 
tudes. The code names of the observatories are given in 
Table 5. This map is based on the nomograph drawn by 


Whaten [1970]. 
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Fig. 50. Magnetograms of the digitized magnetic 
data from 13 observatories. The records have been 


detrended with a 1-5 mHz filter. 
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Table 6 summarizes the polarization character- 
istics of the event at our stations and at the observ- 
atories in the interval 1500-1600. The pulsations are 
Creamy, 10calized in latitudes havingeintensicties 01 
4160y° and 787 0y at SMIT and Churchill respectively. 
(Intensity is defined in Section 2.3.) This localization 
is allustrated in Figure 51. The powers in all three 
components show sharp peaks at SMIT dropping to half their 
peak values only 2-3° to the north or south. At Baker 
Lake (Figure 50) the pulsations have already lost their 
Sinusoidal appearance and at Resolute they are barely 
perceptible. Similarly, in the south at MCMU (Figure 43) 
the pulsations look irregular and at LEDU show clearly 
for only two or three cycles. 

Figure 52 shows the power spectra for our stations 
and the two observatories, Baker Lake and Churchill, which 
are nearest to the centre of activity of the micropulsation 
event. With the exceptions of CAMB and Baker Lake, all 
the H and D component power spectra show distinct peaks 
at 2.1-2.2 mHz. Normally the power in the Z component is 
the smallest but it does reach appreciable values at 
Churchill. Micropulsations recorded at this observatory 


may be influenced by coastal induction since the observatory 
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52. Power spectra of the 3 magnetic components 


dtwourescati1ons, and at baker Lake andeChurchishl. | the 
Spectra were calculated from data recorded in the interval 


1500-1600 UT (day 186). 
ae COMM Ze Twice i. 


The records were detrended with 
The power is in units of y?2 sec. 


(continued on page 129) 
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is situated near Hudson Bay. An indication of a second 
large spectral peak is evident near 3.0 mHz at Churchill, 
SMIT, and FTCH. At CAMB and Baker Lake at least three 
spectral peaks are evident. These peaks are at v 1.5 mHz, 
2.1 mHz, and 3.0 mHz with the low frequency peak most 
dominant at both CAMB and Baker Lake. Micropulsations 

at CAMB have the most complex spectra and the least power, 
with two distinct peaks in the D component and one broad 
low frequency peak in the H component. Based on the 
spectra from our stations and these two observatories, 
there appears to be a dividing line separating a southern 
region with sinusoidal pulsation trains from a northern 
region with complex wave forms and spectra. This line 
must be somewhere near 72°N with Baker Lake directly to 
thesnonthsofutheslinesand SUlteandséhurchilh direct lyato 
the south. As we will see later, this line coincides with 
the demarcation line. 

The latitude-dependent characteristics of the 
powers ,phases, and polarization parameters, associated with 
the 2 mHz Pc 5's are summarized in Figure 51. Data from 
the observatories at Newport and Baker Lake have been 
included-inpthe, plots since thesesobservatories are situated 


near 302° geomagnetic longitude. The powers in all three 
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components reach a sharp peak at SMIT (Figure 51A), 

both the H and D components having comparable powers 

at this station. To the north and south of SMIT the 

D component power decreases the most rapidly of the 

tireess the ellipticities (Figure 510) also show very 

distinGtuslaeutudinalacharacteristics a mine senses. 0. 

polarization in all three planes show distinct reversals 

occurring between SMIT and Baker Lake. From SMIT south- 

ward the micropulsations are polarized in a CC sense in 

the H-D plane, agreeing with the statistical results 

presented Mn aSection 3.2. “Most of stheys tations, have 

CW polarization in the H-Z plane. In the D-Z plane, 

the polarizations are the same as those in the H-D 

plane with stations from SMIT southward having CC polar- 

ization and Baker Lake and CAMB having CW polarization. 

The changes in the senses of polarization in the H-D 

and H-Z planes are apparently caused by an abrupt 160° 

phase change in the H component;,occurring between CAMB 

anadesmun (rigure s518)2 © Thissphaseichangessmalsomeyident 

in the magnetograms of the micropulsations (Figure 43). 
The worldwide distributions of the polarizations 

in the H-D plane are summarized in Bigunvemoct nr Novesthat 


in this figure the sizes of the polarization ellipses 
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Fig. 53. A summary of the latitudinal and longitudinal 
H-D polarization characteristics of the 2 mHz micropul- 
Sations (day 186) plotted in local geomagnetic time. The 
times were determined by using the nomograph drawn by 
Whalen [1970]. The centre of each ellipse is at the 
position of the corresponding observatory or station. 
Ellipses with dotted lines indicate CW polarization and 
ellipses with solid lines indicate CC polarization. 
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are not scaled to the intensity of the pulsations at the 
stations. Comparison of Figure 39 with Figure 53 shows 

that the spatial distribution of the polarization parameters 
of this event agrees remarkably well with the statistical, 
diurnal variations obtained at our line of stations. The 
polarizations at the observatories on coastlines must be 
considered carefully, however, because of possible ocean 
induction effects (Appendix A7). There is certainly 
complete agreement between the statistical diurnal variations 
in the senses of polarization (Figure 39) and the data 

shown here (Figure 53). All the stations and observatories 
south of 72°N and west of 1300 LGT have CC DOTA Zatton . 

Both Baker Lake and CAMB have CW polarization as expected 
and@ierrvodur., which Ws “ast “of 1300 ebGh shaseratner 
elongated but CW polarization. The polarization angles 

at all stations and observatories, with the exceptions of 
thescoas tal opservatories ‘at Sitka. Victoria.,eandms uc. JON. 
are also consistent with those expected from the statistical, 


diurnal variations at our stations. 


(d) Day 195 2200-2230 UT 
This interval shows a preponderance of very high 
frequency Pc 4 micropulsations. Two spectral peaks are 
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at 20.3 mHz. The 20.3 mHz micropulsations have intensities less 
than By" at CAMB, SMIT, FITCH, and MCMU and consequently the 
polarization parameters for these stations may be in consider- 
able error. Figures 54A and 55A indicate that the H component- 
powerseot both the 12.7 and) 20.3 mz micropulsat tonsetrave sharp 
PpeakS=at MENK® The 12.7 mHz pulsations also have ay marked in- 
crease in power in the D component at CALG and in the power in 
Ee we component at CAMB. The latitudinal characteristics of 

ure ellipticities of these two micropulsation trains sare com- 
plicated and do not show any obvious latitudinal reversals in 
Sense. ine micropulsations exhibit ycn polarization inethe Hp 
Dianewatmenessouthern Stations MCMUS MENKSEUEDU ancmGALG. seach 
micropulsation event has two abrupt latitudinal phase changes 


in the H component (Figure 548, 55B). These changes occur at 


2 


65°N and at 70-75°N. The phase changes between SMIT and 
CAMB may not be meaningful since both events have low amplitudes 


and complex spectra at these latitudes. 
(e) Day 262 1900-2000 UT 


This wintervalsas a complicated but interes Ginga one. 
The activity here differs from most other micropulsation 
events in that the dominant frequency in the H component 
decreases Witieincreasing latitude?’™ inere are at least four 
distinct frequency peaks in the H components, and each frequency 
Tomo inant neaananrrow lati cudinaue rancecmmeliiSmued ules is sev — 


dent even from a quick, visual inspection of the magnetograms 
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(Figure 45). This marked latitude-dependence of the H compo- 
nent frequency makes latitudinal plots of polarization para- 
meters based on a spectral peak at one frequency rather meaning- 
less. Consequently we will show no plots of the latitudinal 
polarization parameters, with the exception of those for a 
Os -aisotHz mMicropulsation trainmwntcn as evident at alilane 
Stations. |The power spectra from the 7 stations shows that 
there are 4 distinct spectral peaks in the H components 
(Figure: 56).) A 1.3 maz peak predominates at CAMS (this veak 
does not show in the figures since the data were detrended 
using the 2.5-100 mHz filter), a 4.1 mHz peak predominates 
at SMIT and FITCH, a 6.6 mHz peak predominates at MCMU, and 
a 14.7 mHz peak predominates at MENK and LEDU. CALG does 
not have any distinct spectral peaks in the H component. These 
Spectral peaks and the latitudinal ranges im which they are 
most. dominant are plotted in Figure 57. There is obviously 
a distinct correlation in the latitude-dependence of these 
Pedrseanc tneeStatistical results derivediiarom the intens wey 
Dedks «higure 22). hts observation ts discussed in more 
detat Maite soect1 one 4zale. 

The Z components at MCMU, FTCH, SMIT, and CAMB 
have many of the characteristics of the HK ccmponents, 
showing the same pronounced peaks anc the same latitude- 


frequency dependence. In fact, at FTCH the 4.1 mHz power 
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Fig. 56. Power spectra of the 3 magnetic components 
at the stations. The data are from the interval 1930- 
2000 UT (day 262) and the records were detrended with a 
2.5 -100 mHz filter. The power is in units of y? sec. 


(continued on page 139). 
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Fig. 57., Latitudinal ranges of the dominantepeakss an 
the power spectra of the H components (Figure 56). The 
solid line shows the approximate trend of the intensity 


peaks plotted in Figure 22. 
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in the Z component far exceeds that in the H component. 
At MENK, LEDU and CALG however, the powers in the Z 
components are negligible. The ene enenectre of the D 
components differ completely from the other spectra. A 
Single distinct peakeshows at 9.6 mHz mat CALG. LEDU, 

and MCMU. The behavior of the powers in the D components 
at the northern stations is somewhat more complex and 

no dominant spectral peaks show at FICH and SMIT- 

As already mentioned the 1.3 mHz spectral peak 
occursmat all stations. “Latitudinaleplotsmof. the oowers 
of the three components and the polarization parameters 
associated with this peak are shown in Figure 58. The 
powers in all three components increase monotonically 
with latitude to a maximum at CAMB (Figure 58A). The 
sense of polarization in the H-D plane (Figure 58C) is 
CiM@atecAMB and CC at all other stationsz. In° the H=4and 
H-D planes the senses of polarization are CW at all 
Stations including CAMB. The apparent change in the 
sense of polarization in the H-D plane is caused by a 
gradual latitudinal change in the phase of the D component. 
The phase of the H component changes very little with 
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4. CONCLUSIONS 
4.1 Discussion of the observations 


The number and complexity of the parameters affect- 
ing geomagnetic micropulsations make it almost impossible 
to deduce a simple model which explains all the observed 
features of Pc 4's and Pc 5's. We can only hope to elim- 
inate those theories which obviously conflict with the 
experimental data and perhaps select a few which show good 
qualitative agreement. The most promising source of energy 
for both Pc 4's and Pc 5's appears to be Kelvin-Helmholtz 
instabilities at the dawn and dusk sections of the magneto- 


~_— 


~ Shiow that aa ~ so L-= 
6 SnOw chav Some O1 


pauSe. Evidence exists Citls Sener uy 

is coupled into resonant oscillations of magnetic shells, 
although the simple toroidal mode discussed in Section 1.3 
does not explain the apparent polarizations of the micro- 
pulsations. Tne existence of discrete Pc e4eand Pc.5 peaks 

in the geomagnetic spectrum is probably due to the prefer- 
ential coupling of micropulsation energy into magnetic 

shells near the auroral oval and near the plasmapause. At 
present the spatial characteristics of the micropulsation 
fields are too poorly understood to explain with any accuracy 


the observed polarizations in the H-Z and D-Z planes. 
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Ine the discussion of micropulsation intensities in 
Sect lOMmesne 1 We noted that the intensity maxima of Poe 51s) to le 
HOW thes euroral oval (Figure 19). “This observation is claritied 
Pie ture soo whic SHOWS ad polar plot Of thesnOSa1 clonseon ste 
intensity maxima, together with the auroral oval and the aver- 
age position of the plasmapause [Carpenter, 1966]. Except 
around local noon, the auroral oval and the zone of. maximum in- 
tensity overlap. At local noon, however, the centre of activity 
has obviously moved north of 70°N and unfortunately we have no 
stations between 68°N and 77°N. Consequently the mismatch in 
the positions of the auroral oval and the activity centre may 
not be real. Another interesting feature is the proximity of 
the plasmapause to the auroral oval in the interval 2000-2200 LGT. 
This bulge in thé plasmapause could considerably affect the pro- 
pagation of the dayside micropulsation energy to the nightside 
of the earth. The bulge would certainly complicate any latitude- 
dependent characteristics of micropulsations occurring on the 
duskside of the earth. 

Fairfield [1968] has shown that the dayside segment 
of the auroral oval maps along geomagnetic field lines to the 
magnetopause. In addition, Frank [1S AlMebasechownmchaterieald 
finessatrethese frghelatitudes#iierveryecioce tothe daysude 
polar cusp. Both these features indicate that the energy 


of low frequency Pc 5's might be derived from 
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Hid 259. ePOlarm plots of themwauroraimove lm ibeLween 
dotted lines), the intensity maxima of Pc 5's (between 
solid yes and the plasmapause (dashed line). The 
auroral oval corresponds to the oval deduced b } 
[1967 | for periods with moderate activity (Q ee ete 
circular grid lines show the corrected geomagnetic 
coordinates. The time scale is marked in hours of correc- 


ted geomagnetic local time. 
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an interaction of the solar wind with the earth's magneto- 
sphere. We will show later that our polarization studies 
Suggest that this interaction sets up a Kelvin-Helmholtz 
instability near the magnetopause, and that this energy 
1smguided along field lines. to the earth’ sesurface. self 

the energy does indeed come from the solar wind, we still 
have the problem of explaining why the centre of activity 
follows the auroral oval on the evening side of the earth. 
One possible explanation is that the micropulsation energy 
is coupled into the “inner dipole" region of the geomagnetic 
field and that this whole region oscillates in response to 
the energy input on the dayside. Since the auroral oval 
lies near the region of the transition from open to closed 
field lines and the beginning of the plasma sheet [Vasyliunas, 
1968], any micropulsation energy north of the auroral oval 
on the nightside would not be expected. This is especially 
true since the plasma sheet effectively absorbs micropul- 
sation energy [Barnes, 1966; fee, Ware, 

An alternative explanation for the micropulsation 
energy in the nightside auroral oval is that the energy 
originates in the magnetotail and reaches the earth by 
propagating along field lines [Rostoker, 1967]. Tamao [1968] 
and Saito and Sakurat [1970] suggested that this energy 
might be derived from hm instabilities near the inner 


boundary of the plasma sheet. 
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In the plots of intensity versus latitude (Figure 20) 
we noted that Pc 4's have proportionally greater intensities 
at low latitudes than Pc 5's. The lowest frequency Pc 5's 
(0-2 mHz) have maximum intensities near 75°N whereas Pc 4's 
(10-20 mHz) have maxima near 60°N. This frequency -dependent 
transfer of energy to lower latitudes suggests that at least 
some of the energy of the micropulsations is distributed in 
resonant oscillations of magnetic shells similar to the 
toroidal mode discussed in Section 1.3. Micropulsations 
at CAMB have considerable energy in all spectral bands, 
indicating once again that the field lines in this region 
are near the source of micropulsation energy. 

The latitude-frequency dependence of the intensity 
peaks is shown better in the scatter plots of Figures 21 
and 22. Both these figures show distinct latitude-frequency 
trends in the groupings of the intensity peaks. Figure 22 
also shows a substantial cluster of peaks covering all 
frequencies from 1-20 mHz at latitudes above pe RS 
high latitude cluster indicates once again that stations 
at these latitudes lie near the source of dayside micro- 
pulsation energy. This feature does not appear on the 
nightside plot (Figure 21) and most of the intensity peaks 
aré south of 72°. The four Pc 4 peaks at CAMB (77°N) may 
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simply indicate energy input from the tail region and may 
not have any connection with micropulsations at lower 
latitudes. 

The second feature evident on both plots is the 
latitude-frequency correlation of the intensity peaks 
south of 75°. This feature certainly indicates that 
micropulsation energy is distributed in resonant oscil- 
lations of magnetic shells. The frequency estimates in 
Figures 21 and 22 are higher than the fundamental toroidal 
eigenfrequencies calculated by Dungey [1954] for a normal 


dipole field in a constant density plasma ( pe = i oe cm73) 


but agree better with the eigenfrequencies of a compressed 
dipole field in a cavity of radius 8 Re with a plasma 
HSE Ole OAs in) 2" g cm-> [Westphal and Jacobs, 1962]. 
At this point we might ask which component, if 
any, shows the most distinct and predictable latitudinal 
dependence in its intensity peaks. Actually the particular 
component seems to depend both on local time and the frequency 
OpeLReTGIVeNn MmiCrOpULSdtiONs. | BOUNsuncatmdanG  UecOmponentus 
Ofmec oes Show strong iatitudinal peaksein the morning. 
but normally the H component shows the most distinct 
latitudinal peaks in the afternoon (see for example 
Pigunesmao dndeol)., re 47s on the other hand often show 


strong latitudinal peaks in the D component in the morning 


(see for example Figure 46). 
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These features certainly do not describe a 
SiMplemtOGOnad! OSCiuliat1 On wei adcitioneeties limited 
longitudinal extent of most pulsations precludes the use 
of the toroidal equation for a description of the modes 
of oscillation. The asymmetric,guided poloidal mode is 
also not a good approximation since it predicts linear 
polarization in the H direction, and both Pe 4's and 
Pc 5's normally have elliptical polarizations in the 
H-D plane. 

We might intuitively expect that since the 
Symmetric toroidal mode is polarized in the D direction 
and the highly asymmetric poloidal mode is polarized in 
the H direction (at the earth's surface), more generalized 
modes of resonant oscillations of magnetic shells would 
be elliptically polarized in a plane perpendicular to the 
main field. Cummings et al. [1969] have reported micro- 
pulsations with this type of polarization in magnetic data 
FrOtethessadtelliteeAlso i.) TN1S Sate lincemis winwaneequa lon 
eames V:Di card Lao 200 Re Situated on field lines which intercept 
the earth's surface at v 67°N. The micropulsations detected 
by the satellite occupied two frequency ranges, one centred 
at 5.3 mHz and the other at 9.8 mHz. The frequencies of 
these spectral bands agree qualitatively with the frequen- 
cies of the micropulsations we find with intensity peaks 


at 67°N. 
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Since the micropulsations observed at ATS 1 
were transversly polarized we should perhaps consider 


the hm wave equations with the additional constraint 


where by is the perturbation magnetic field in the direc- 
tion of the main field. Cummings et al. [1969] and 
Wamgaladze and Brunellt [1970] have considered in some 
detail the resulting simplifications in the wave equation 
resulting fromatnis constraints @fhesconstrain by = 0 
(using the coordinate system introduced in Section 1.4) 


and Faraday's Law 


indicate that 
) ) y 
Bx (h3Ea) a eae sees 0 (eile) 


Consequently equations (3) and (4) become 
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Equation (12) is equivalent to the equation for the 
toroidal mode: equation (13) describes the highly asym- 
metric or guided poloidal mode. Both these equations 
describe the resonant oscillations of magnetic shells 

and there are derivatives only along the field lines 
(ive. ws, To solve these wave equations we must 

once again apply the boundary condition E = 0 at the 
earth's surface. We then find that the modes can exist 
only for certain discrete frequencies, and that in both 
cases the fundamental frequency decreases with increasing 
latitude. Cumming et al. determined, however, that the 
fundamental frequencies of Eo and E3 are different. 
Thus equation (11) can never be satisfied for the funda- 
mental mode. This they construed as meaning that the 
fundamental mode cannot exist. They found, however, that 
all the higher harmonics of E. and E3 have almost 
equal frequencies. In other words, if the harmonics of 


2 and E3 are w. and 0, respectively then 
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If the plasma density in the plasma trough varies 


as r and if the equatorial number densities are in the 


range 1-5 elfen. then according tomtnescalculations: oF 
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Cummings et al. the expected frequencies of the second har- 
monics occurring at 67° are in the range 1/-38 mHz. These 
frequencies are much higher than we would expect from our 
datas For the frequencies to match those in our data the 
equatorial densities of electrons would have to be 10 - 30 
el/em>, Mmeserdensities are clearly mucns too m10ne sin 
fact, the calculated frequencies of the fundamental modes 
compare best with our data. 

Up to the present there has been considerable 
confusion regarding the possibility of a latitudinal de- 
pendence in the frequencies of micropulsations. Our analysis 
Shows that with few exceptions a micropulsation train has the 
same frequency at all latitudes, although the latitude of 
the peak intensity tends to decrease with increasing fre- 
quency. The earlier work of Fllis [1960] and Obertzg and Ras- 
popov [1968] confirm these observations. Conversely, many 
Reseancners have found. through statistical anatysis, that the 
average frequency of micropulsations decreases with increasing 
Hatitude.. (Compare Figure 3 with our datesin Figure 27.) 
Phesesstatistical results dOenot NECeSsani | Ve contrada cusOUr 
observations since those events which have a frequency equi- 
valent to the resonant frequency of a given magnetic shell 


would typically have thier greatest amplitudes at the latitude 
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of that shell. Averaged over many events the larger 
amplitudes at this resonant frequency would tend to give 
more weight to events with this frequency. Since the 
resonant frequencies of the magnetic shells decrease with 
increasing latitude, a latitude-dependence in the average 
frequencies of the micropulsations would be evident. 

In special cases, when considerable broad band 
energy is introduced into the magnetosphere (e.g. a sudden 
impulse), there might be a pronounced latitudinal depen- 
dence in the frequency of a given micropulsation event. 
For example, suppose the input energy has spectral peaks 
at 3 distinct frequencies fi, fo, and f, (Figure 60). This 
energy is coupled into resonant oscillations of magnetic 
Shells with energy of frequency f4 at the lowest latitude 
and fy at the highest. Consequently micropulsations in 
region c (Figure 61) have an apparent frequency Fy. in 
region b_ they have an apparent frequency Fo > and in 
region a an apparent frequency f3. Thus there is an 
evident decrease in frequency with increasing latitude. 

In reality the situation is more complex since 
the spectral peaks of the input energy are controlled not 
only by the source mechanism, but also by the size and 
shape of the magnetosphere. The boundary at the outside 


of the magnetosphere forces the whole magnetosphere to 
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Fig. 60. An example power spectrum of a broad band 
Source of micropulsation energy. 
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Fig. e619) alhe latitudinal: distnibuti onmof = thesenergy 
and frequency of the micropulsations. 
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oscillate at preferred frequencies. In an axisymmetric 
magnetosphere the modes of oscillation would be very 
similar to the simple poloidal mode [Kato and Akasofu, 
1955; Jacobs and Westphal, 1964; Carovillano et al., 
1966]. This poloidal energy is coupled into resonant 
magnetic shell oscillations with the lowest frequencies 
near tne auroral oval? 

The event occurring on day 262 shows a distinct 
latitudinal dependence in the frequency of the H_ com- 
ponent and may, according to the above discussion, be 
caused by a broad band energy source. This event is 
similar to the micropulsations for which Voelker [1968] 
found latitude-dependent frequencies in the H components. 
In fact the latitude-frequency dependence found by Voelker 
compares remarkably well with our data (Figure 57). Since 
these events have a pronounced latitude-frequency depen- 
dence only in the H component, they could be manifestations 
of the guided poloidal mode. Equation (13) describes the 
perturbation electric fieldstor this modeseeCummingseet als, 
in determining the resonant frequencies of this mode, have 
shown that the harmonics have the same frequencies as the 
harmonics of the toroidal mode. Even the fundamental 
frequency of the guided poloidal mode is not too different 


from that of the toroidal mode (e.g. 17 and 12 mHz respec- 
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tively at 67°N for a rot distribution of the density of 
electrons). Consequently we would expect the guided 
poloidal mode to show essentially the same latitude- 
frequency dependence as the toroidal mode. 

Stebert [1964] proposed another explanation for 
the latitude-dependent frequencies in the H components 
of some micropulsations. He postulated that latitudinal 
changes in frequencies are a result of field-aligned hm 
oscillations in the magnetosphere. For these oscilla- 
tions to occur, the magnetospheric plasma density must 
have a lamellar structure, with the density changing 
continuously along geomagnetic field lines and changing 
abruptly in directions perpendicular to the field lines. 

In any case, the energy for the latitude-dependent 
oscillations probably comes from a sudden compression of 
the dayside magnetosphere by the solar wind. Voelker has 
shown that SSC's and Si's associated with pulsations with 
latitude-dependent frequencies occur near local noon (as 
does our event), thus giving additional support to the 
compressional model. 

If the energy of Pce4candshceoemicropulsations 
originates on field lines in the auroral oval and is coupled 
into resonant magnetic shell oscillations at lower latitudes, 


it is difficult to see why there should be distinct Pc 4 
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and Pc 5 classifications. Satto [1964] and Campbell [1967] 
have shown, however, that there is a natural division of 
the average geomagnetic spectrum into Pc 4 and Pc 5 bands. 
Both these authors have computed spectra which show sub- 
SboneualeneakS ine Dota tie pCu4mandseceompands ame ln med dc ita Onn 
inspection of Figure 22 shows that in our data there is an 
obvious concentration of micropulsations in the Pc 4 band, 
most of the intensity peaks being south of 65°, 

The division of the geomagnetic spectrum into 
Pc 4's and Pc 5's possibly indicates that these micropul- 
sations have different sources of energy. Pc 4's, for 
example, might derive their energy from a bounce resonance 
interaction of particles and waves (Section 1.3), whereas 
most Pc 5's derive their energy from instabilities at the 
magnetopause. The sinusoidal appearance and long wave 
trains of many Pc 4's make the bounce resonance tneory 
particularly appealing [Dungey and Southwood, 1970]. On 
the other hand Pc 4's and Pc 5's often obey the same diurnal 
polarization rules (see for example Figure 24 and Rankin 
and Kurtz [1970])suggesting that they have the same energy 
sources. 

AheelormawLonsol GisGincesy Gedmaldm) Geoms DOC.ciad | 
peaks might also result from the preferential coupling of 


energy into given magnetic shells, say at the plasmapause 
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and near the auroral oval. Hirasawa and Nagata [1966] 
proposed that Pc 4's are caused by a strong coupling of 
energy into resonant oscillations of magnetic shells near 
the plasmapause. Tamao [1968, 1969] suggested that a 
trough in the spatial distribution of the Alfvén velocity 
serves aS a trapping region for hm waves. Such a trough 
exists just inside the plasmapause [Dungey, 1967] and we 
might expect a large number of micropulsations with 
frequencies equivalent to the resonant magnetic shell 
frequencies in this region (~ 10-20 mHz). 

Figure 62 summarizes the diurnal polarization 
characteristics of Pc 5's in the H-D plane and should be 
compared with Figure 39. Figure 62 shows the polarization 
characteristics of micropulsations with mean frequencies 
near 5 mHz. The features of the diagram must be moved 
to the south for micropulsations with higher frequencies 
and to the north for lower frequencies. Actually the 
diagram is probably a good indicator of the characteristics 
Of polarizations Of PCce4)syasmwelle | ines verticals linesat 
2000-2100 (~ 1130-1230 LGT) marks the temporal reversal 
noted simultaneously at all our stations [Samson et al., 
1971]. This reversal also shows in the data of Kato and 
Utsumt [1964]. A second reversal at approximately 0300- 


O400m (1830-1930 UG) is at thesedgemotethesdiagram: 
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The short period (7-8 hr) of CW polarization occur- 
Gingeat tne "southern stations in tite sarternoon sugdests sthat 
ciewdayeve.divrded IntG 3 sectors face mariing (0700-1200 LG). 
afternoon (1200-1900 LGT) and night (1900-0700 LGT). The night- 
side and morning sectors have the same polarization configura- 
tions but different energy outputs. Impulsive Pi micropulsa- 
tions dominate much of the activity in the nightside whereas 
continuous Pc micropulsations dominate most of the activity in 
the morning. 

A plausible explanation of the polarization reversal 
dusilsG-1730 Gh "rs that the energy of “daytime micropulsdtions 
originates from a Kelvin-Helmholtz instability at the magneto- 
pause [Atkinson and Watanabe, 1966; see also Section 1.3]. The 
Sensei polarization OF Pc S's "occurring ates vat OnswsOouUth OF 
CAMB is CC in the morning and CW in the afternoon in agreement 
with the sense predicted by this theory. The Varge tnput of 
energy at all frequencies on the dayside auroral oval (see Fig- 
Uveeweey oives further Support to this vneory = Unrortunatery 
this model does not explain the latitudinal change in the sense 
Ofap0larazaciron. 

The demarcation line lies in the curved region above 
the line of maximum intensity. (The demarcation line has been 
labeled demarcation region in Figure 62 because of the inaccur- 
Acyeinedetermimning the position of this @line with respect “torthe 
latitude of the maximum intensity.) The second line (dotted 


curve) at high latitudes is necessary Since certain 
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times exist during the evening when Pc 5's at all our 
Stations, including CAMB, have the same sense of polar- 
ization in the H-D plane. Unfortunately the two northern- 
most stations are 9° apart and the region of reversed 
polarization (directly above the demarcation line) is 
sometimes not evident in the evening data. 

At first we may be tempted to associate the 
polarization demarcation line with the position of the 
auroral oval, especially since most low frequency Pc 5's 
have their maximum intensities and their demarcation lines 
in this region. This feature does not always appear to 
be the case, though, since there is a close correlation 
between the frequency of a micropulsation train and the 
latitude of its demarcation line (Figure 34). Many Pc 4's 
have demarcation lines well south of the auroral oval 
(see for example Figure 46), although in general Pc 4's 
have rather complex latitudinal polarization characteristics. 
Clearly more stations are needed in the region between 68°N 
and 77°N in order to compare the latitudes of the demarcation 
lines of Pc 5's with their intensity centres and the auroral 
oval. 

The station-to-station phase comparisons indicate 
thateabrupes stathtudina lephases changes inwvertnereunealimor 
D component may be more useful than polarization reversals 


in determining the position of the demarcation line. In 
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some cases the relative latitudinal phase changes in these 
components are very large but not large enough to cause a 
change tn the sense of polarization. 

As mentioned in Section 3.2.2 there is a strong 
diurnal variation in the frequencies of micropulsations 
with demarcation lines at a given station (see for example 
Figure 33). Figure 63 is a composite of diurnal polar- 
ization data from SMIT and FTCH and shows the sense of 
polarization as a function of time and frequency. The 
dotted line indicates the frequencies of micropulsations 
with demarcation lines near these stations (we shall call 
this the demarcation frequency). This frequency is as low 
as 4 mHz in the interval 1000-1200 rising to a maximum of 
14 mHz from 1800 to 0400. Evidently the diurnal variation 
in the demarcation frequency is closely associated with 
theslatitudina Wmottonso® the jauronal #ovaleeeeAteoures tations 
the oval is furthest south between 0800-1200 (2330-0330 LGT) 
and moves quickly northward at approximately 1500-1700 (0630- 
0830 LGT) (see Figure 59). Although not evident in our data 
a diurnal trend in the demarcation frequency is apparent in 
micropulsations recorded at LEDU [Kurtz, 1969]. (The 
frequency~dependent reversal is not noticeable in our data 


since our system cannot detect many micropulsations with 
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Fig. 63. The frequency-dependent and time-dependent 
changes of the H-D ellipticities of micropulsations 
recorded at FTCH and SMIT. In the dotted regions, the 
ellipticities are less than -0.1 and in the lined regions, 
the ellipticities are greater than +0.1]. The demarcation 
frequency is indicated by the dotted line. Spaces with 
an N have no data. 
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frequencies greater than 20 mHz). Figure 64 shows a 
comparison of the diurnal variations of the demarcation 
PGcCcuenCHCsmd tat eUUmandeat soil lel CH mmmeVidentlyeetie 
demarcation frequency at LEDU is much higher than that 
at SMIT and FTCH. The difference in the demarcation 
frequencies would be expected if the latitudes of the 
demarcation line were associated with the intensity peaks. 
Since many Pc 4's have intensity peaks near and south of 
LEDU, we could expect that LEDU would show ellipticity 
reversals in the Pc 4 spectral band whereas SMIT and FTCH 
would show reversals in the Pc 5 band. The diurnal varia- 
tion of the demarcation frequency at LEDU also does not 
have the same shape as that at SMIT and FTCH. This might 
Demdmenes Ul PwOtetiesc LOSE, DLOXIMicyeO1 BLEDUMULOBnPelLGmai mes 
following the plasmapause (from 1800-0100 LGT, LEDU is 
actually south of the expected position of the plasmapause). 
The dotted line in Figure 64 shows the high 
frequency limit of the two spectral bands noted by dirasawa 
and Nagata [1966] and Wagata and Fukunisht [1968] in their 
analysis of micropulsation data from Kakioka (36°N) (see 
Sectu1om isc) theyeattriputed thesspectralmband centred 
at 15 mHz to eigenoscillations of magnetic shells near the 
plasmapause since the diurnal variation in the mean frequency 


of the band followed tne diurnal variation in the latitude 
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Fig. 64. The diurnal variation of the demarcation 
frequency at SMIT and FTCH (upper solid curve) and at 
LEDU (lower solid curve). The dashed line shows the 
high frequency limit of the two spectral bands noted 
by Nagata and Fukuntsht [1968] during periods of low 
magnetic activity (K, v 0-1) and the dotted line shows 
the upper limit during periods of moderate activity 
(K, ~ 2-3) : 
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of the plasmapause. On the other hand, they attributed 
the 30 mHz band to hm modes trapped between the plasma- 
pause and the ionosphere. An obvious correlation exists 
between the demarcation frequency at LEDU and the two 
spectral bands noted by the above authors. It appears, 
then, that the ellipticity reversal at LEDU is related to 
the Pc 4 resonance near the plasmapause whereas the ellip- 
ticity reversal at SMIT and FTCH is related to the Pc 5 
resonance in magnetic shells near the auroral oval. 

The polarization demarcation line and the assoc- 
iated frequency-dependent senses of polarization now appear 
to be fundamental characteristics of magnetospheric modes 
OfeeOS Cin dt aiON Odie? GS. and P Ce5.S memicoe cea 0, een 
solving the problem of coupled toroidal and poloidal modes 
of oscillation in a hm wedge, found that resonances appear 
on given magnetic shells and that the sense of polarization 
of the perturbations changes across the resonance region. 
The whole cavity oscillates at a given frequency and both 
toroidal and poloidal modes are evident throughout the 
cavity. The toroidal modes, however, have spatial resonances 
One Givens maghetces hel lSeer hes DOSmtiOnsmomeLnies ems ells 
are determined by the eigenfrequencies appropriate to the 
shells. The polarizations derived from this model agree 
at least qualitatively with the latitude-dependent senses 


of polarization found in our data. 
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The diurnal variations of the polarization 
angles in the H-D plane, shown in Figures 36 and 39, 
indicate that the ellipses of Pc 5 micropulsations are 
oriented in the H direction throughout most of the day. 
It might be conjectured that this is an ionospheric 
Grnecr. = perianps due to =the lati tudinalmlocahizat! Oneor 
the micropulsation field [Inoue and Sehaeffer, 1970]. 
Satellite data indicate, however, that polarizations in 
the H direction are characteristic of the magnetospheric 
modes of propagation. Cummings et al. noted that sinu- 
SOrdalemreropuiSations recorded) at vAls-imnormal ly enad much 
larger radial than azimuthal field components (see Figure 
9). In their data the average orientation of the major 
axis was 30° east of the radial (outward) direction. Patel 
[1965] also noted a preference for greater amplitudes in 
the radial direction near local noon. He found, however, 
that the azimuthal component of the pulsations increased 
toward the dawn side of the magnetosphere. His observations 
agree with our data since we find elliptical polarization 
in the early morning with a preference for polarization 
in the D “direction at SMIT “and FICH in the “interval 1600- 
700010 73021) 130 Gl). "(See 'also Samaontetear. LISTE). ) 


Wisneiiaise 28 te heat aT 
ane a” 

Re ae n@one. % ve n . 

aus anoigertyaeisin & 34 %o aan 8: 


7 4 


ie ae Yo Prd Soodpenidd ROPE We ee a r | 
jie 

a ee eelten ys ‘¥ 

yo nagfexttese! foatoustest ole Oo: bgt eqedtee ot 7 

ebay Ads andoe foc eet) BIEN vobreetuqorstel #9 

nt sgaties eetoc 75a ' 4eveaga ,steothnt evab Ost By i 


sPatitotengsn <i2 10. 9 Tetvalosrens’ w5 “oe 
suats 384) boson fo 29° tenéeney -.o8enesgotd) fa eee W 


Ssuinet yil@mon | CTA te Vaiaepes cweltwataqora te tabtee 

erp Gry) stevowgme> UfePy Dewduwtge asad (oteey saat 
véthe add ty soi cater tio gpehewn od B9eb 4 Tene ee 

Yake® (nePioeeth (rewind) lather ea) te seem SOR eee -_ 


wt Qabud (few oss0e7¢ "OF saupistew @ Seton ote caaer iv 
-edugmod_.beve? 40 .vo0d. Paget seed unt say th Tots wr odd Fi 
Bateevseh anstese! us ott % ifosrqne> tetfomtse ‘3 | 
ena tiavidddh 2th SWdged Sopa 30) ‘to sh?s wed ‘a6t b A WO 5 
notfestunton foatigtl ie wert ob Ssate egab sue a | 7s 
aptaes jteroq nt -oitaiary ¢ Ob tw perotom qin Ms a nt 


7 wet ne oo eee hie 


168 


In addition to being the time of the most circular 
polarizations, the early morning also has the greatest occur- 
Felcemi requency #01 7c aol) Cropu lsat Onsm(seemocGr Tolmom: aces 
Field lines which intersect the earth's surface in the early 
morning sector are swept back from the sun and meet the 
magnetopause at a position where the solar wind flows tan- 
gentially ([Fatrfield, 1968]. Since the magnetopause is 
most susceptible to a Kelvin-Helmholtz instability in this 
region [Southwood, 1968], we would expect most of the 
micropulsation energy and most circular polarizations in 
this region. 

If the energy for daytime Pc 4 and Pc 5 micropul- 
sations originates from a Kelvin-Helmholtz instability at 
the magnetopause, it is difficult to explain the apparent 
difference between micropulsation activity occurring in 
the morning and in the evening. The morning sector usually 
has the greatest number of sinusoidal pulsations and has 
consistent. sometimes almost circular, CO polarization. 

The afternoon sector, on the other hand, usually shows very 
few sinusoidal pulsation trains, a mixture of polarization 
senses, and a predominance of almost linear polarization 
in the H direction. Intuitively we would expect that a 
Kelvin-Helmholtz instability would be equally probable in 


the morning and afternoon sectors of the magnetopause and 
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that this instability would produce micropulsations with 
equally elliptic polarizations but with opposite senses 

in these two regions (see Figure 11). The actual asymmetry 
in micropulsation activity might be caused by the attenu- 
ation of the CC or right handed mode in the ionosphere or 
magnetosphere [Field and Greifinger, 1965; Inoue and Schaeffer, 
1970]. Cummings et al. found, however, that most of their 
events occurred before 1200 LMT, in agreement with ground- 
based observations of the occurrence frequency. Further 
comparisons of satellite and ground-based data should 
Clarify this point. 

One other possible cause of the asymmetry in 
micropulsation characteristics is the dawn-dusk asymmetry 
in the plasma sheet [Vasyliunas, 1968b]. The plasma sheet 
is well defined near the dusk magnetopause whereas in the 
morning sector plasma densities change continuously from 
the plasmatrough to the magnetopause. A sharp discontinuity 
in the density of the plasma might result in considerable 
absorption of micropulsation energy in the evening sector. 

The Z components of micropulsations observed on 
the ground are undoubtedly due to spatial gradients in the 
currents induced in the ionosphere by hm waves. It is 
illuminating to consider current systems compatible with 


the H and D components recorded on the ground and to compare 
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the predicted Z variations and polarizations with those 

we have observed experimentally. Wilson [1966] used a 
rotating line current as a model for the current systems 
induced in the ionosphere by hm waves. In most cases this 
model does not explain our observations, especially the 
variation in the sense of polarization in the H-D plane 
with latitude. This model also does not predict a large 

Z component at the latitude of the maximum in the intensity 
of the horizontal components, which we often find to be 
the case in our data (see for example Figures 47, 50, 53). 
Eleman [1966] also noted this discrepancy between Wilson's 
model and the intensities of the Z component. 

The large latitudinal phase shift occurring in 
the two horizontal components is indicative of a rather 
strong shear flow in ionospheric currents. Since the phase 
shift is usually most noticeable in the H component [Kaneda 
et al. 1964; Obertz and Raspopov, 1968], especially during 
morning Pc 5 activity (see for example Figure 51), the shear 
must occur during that part of the pulsation cycle with 
Strongseast-west currents. = [hisesnear iseeyidenteingthe 
equivalent current system for Pc 5's given by Jacobs and 
Sinno [1960]. Normally, however, the current patterns must 
change during the pulsation cycle in order’to produce the 
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Consider now a micropulsation event with nearly 
circular polarization in the H-D plane. The latitudinal 
dependence of the phases and powers are shown in Figure 
65. (Compare Figure 65 with Figures 51B, 55B.) Since 
there is a very rapid latitudinal phase change in the H 
component we would expect an equivalent current system 
Similar to that shown in Figures 66 and 67. The moving 
current vortices proposed by Obertz and Raspopov could 
possibly cause the changes in direction of the current 
flow shown in Figure 66. However, it is not easy to see 
how their current vortices could produce sinusoidal oscil- 
lations. We will call the region of shear flow the 
demarcation line because of its association with the polar- 
ization reversal. Directly to the south of the demarcation 
line the Z component is 180° out of phase with the H con- 
ponent (Figure 67) and directly to the north of the demarcation 
line the Z component is in phase with the H component. The 
maximum amplitude of the Z component occurs very close to 
the demarcation line. Thus we expect linear polarization in 
the H-Z plane with the value of the polarization angle depend- 


ing on whether the station is to the north or south of the 


demarcation line. 
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Rig-s OS a latitudinal plots=ot thespowerssandsphases 
of the H and D components of a model micropulsation 
event. 


lu =. o° i i ae 
Q ——> 
ST een it pac rmaaae, iL tb 
< lies He 0, Ge 
deme 


Fig. 66. The equivalent ionospheric currents of the 
model micropulsation event (viewed from above). This 
diagram shows the equivalent currents during two parts 
Gf ethe «cycle wii chi@are, 90° Papart. Siiieedottedecincies 
showethe’ direction of ithe rotatvonweotf the tcurrentavectors:. 
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GEOMAGNETIC LATITUDE 


Pid..0/. the magneticuperturbationss Of sthenueand. 7 
components of equivalent currents flowing to the west 
(dotted region) and to the east (lined region) in the 
ionosphere. This figure corresponds to the 0° equival- 
ent current system shown in Figure 66. 
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Further south of the demarcation line, the Z 
component is again in phase with the H component (region 
a Figure 67). In this region the polarizations in the 
D-Z plane should be opposite to those in the H-D plane. 
This is evident in Figure 30. All the stations which are 
Faresouth of the centre of activitystomerc os MCMUss MENKs 
LEDU, and CALG) have polarizations in the D-Z plane which 
are essentially opposite to those in the H-D plane. This 
observation is particularly evident if we consider the 
diurnal, trends ain the ellipticiticssesihe ellipticities 
in the H-D plane have a U-shaped diurnal variation, whereas 
in the D-Z plane they have an inverted U-shaped variation. 

The above model is of course an over simplification. 
In general the polarizations will be highly elliptical and 
the actual demarcation line will run at an angle to the D 
direction. We also have not considered three-dimensional 
current models such as those proposed for polar magnetic 
substorms [Bonnevter et al. 1970]. Our model does suffice, 
however, to explain some of the more evident latitude - 
dependent characteristics of Pc 5's and, to some degree, 
those of Pc 4's. In order to establish this current model 


for Pc 5's we need more magnetic stations between 68°N and 
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The analysis of the model of ionosphere currents is 
considerably simplified by computing instantaneous, lati- 
tudinal profiles of the amplitudes of the three components. 
Inese can be calculated directly from the cross spectral 
estimates Pg shea We have computed amplitude profiles 
for a few of the events in this thesis but there are usually 
too few stations north of the demarcation line of Pc 5 events 
to clearly establish the latitudinal dependence of the com- 
ponents. The latitudinal profiles of Pe 4%s are generally 
very complex possibly because the Z components are strongly 


affected by induction. 


4.2 Summary 


From the analysis of the data in our experiment, 
several features of Pc 4's and Pc 5's have been determined. 
First, most micropulsation trains have frequencies which are 
independent of latitude. However, micropulsation trains with 
apparent latitude-dependent frequencies in the H component do 
occur. The guided poloidal mode appears to satisfactorily 
describe these latitude-dependent micropulsations. 

ht Meast some of the energy v0terce4 and Pc Sonicro- 
pulsations is distributed in resonant oscillations of mag- 
netic shells. However, the simple toroidal mode does not 


adequately describe the spatial features of the micro- 
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pulsations. The latitude-dependent frequencies of resonant 

oscillations of magnetic shells are evident in two ways. 
(a) The latitudes of the intensity maxima of 
Spectral components of micropulsations decrease 
with increasing frequency. 
(b) The sense of polarization of micropulsations 
changes across the demarcation line. The 
latitude of the demarcation line decreases as 
the frequency of the micropulsations increases. 
Actually an abrupt latitudinal phase shift in 
either the H or D component is probably more 
useful than the polarization change for deter- 


mining the latitude of the demarcation line. 


Thesexistence of distinct, PG 47and) PC. 5s peakse in 
the average geomagnetic spectrum is probably due to the 
preferential coupling of energy into magnetic shells near 
the auroral oval and near the plasmapause. Pc 5's occur 
predominantly in the auroral zones whereas Pc 4's occur 
predominantly near the plasmapause. 

The reversal in the sense of polarization which 
occurs between 1130-1230 LGT in the H-D plane gives strong 
evidence in support of a Kelvin-Helmholtz instability as 


the source of energy for daytime Pc 4's and Pc 5's. The 
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argeminimlux@ofepowersateal | frequencies in the daytime 
auroral oval further supports this theory. Unfortunately, 
the theory does not explain the apparent differences 
between the morning and afternoon micropulsation activity. 
It is also difficult to see how any hm instability could 
produce thesvery* sinusoidal’ Pcedts which» often occurein 
ouneda tas 

The Z components and the related polarizations 
in the H-Z and D-Z planes of micropulsations observed at 
the ground are undoubtedly caused by spatial gradients. 
in the ionospheric currents associated with the micropul- 
Savionss “During part of. the pulsattonwecycie: sstrongasnear 
currents flow to the east and the west in the ionosphere. 
Unfortunately, the wide spacing of our stations in the 
northern regions and the complexities in the diurnal 
variations of the polarization parameters have limited 
our interpretation to a very simple ionospheric current 
model. An accurate current model would be very useful 
in relating ground-based observations to the hm modes in 
the magnetosphere. 

There are a number of obvious ways to supplement 
the data collected in this experiment. In the region 
between 68°N and 77°N a closer Spacing of stations would 


help in comparing the relative positions of the intensity 
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maxima and the demarcation lines of Pc 5's. A closer 
Spacing would also help to develop a more accurate model 

Gia thesvonospheric current sys temseot PCa uSeu eo incesour 
fluxgate magnetometers are inherently limited to micro- 
pulsation frequencies below 20-30 mHz, it would be helpful 
to complement our fluxgate records with data obtained from 
OLNC re inSit RUMEN GSS UC -dSmiNGUctl ON mcOlls a ahi Nnadeny. 
stations positioned along an east-west line are necessary 
in order to compare the phases of the H, D, and Z components 
at different longitudes, and to compare longitudinal polar- 
ization characteristics with the diurnal characteristics, 
found in this experiment. These stations must be established 
in regions with uniform geoelectric characteristics, other- 
Wise GCarthvanduction effects could mask the atruess ounce 
configuration of the micropulsations. Knowledge of longi- 
tudinal phase changes would allow the development of a more 


accurate model of ionospheric currents. 
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ABP END Uke Ale ANALY SISMOF COMPEE AS VECTORS TIME SERIES 


Al.1 The general correlation matrix 


A time dependent vector field is completely 

specified by the set of time series ei stile x5 (7%, t), 
x,(¥, t) where X1> Xo» X32 are components of the vector 
field measured in 3-space, y the Spatial position of the 
measurement extending over the whole vector field, and 

t the time covering the duration of the event. In general 
these series are difficult to interpret and we find it 
useful to introduce a different representation of the data 
in the form of a correlation matrix. In three dimensions 
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(Al) 
where * denotes the complex conjugate, r and iby are varying 
spatial positions, and t the lag [Wolf, 1954]. The expecta- 
tion ~ 


lim 1 a 
Too 27 x; (t)x5(t-t) be ACN) 


ai 


siete). x5 (t-1)> - 


wave 


er fein 


ifs2 7 r ¢ ' 
r! b- here eh Caos “ss sanen 
= 2 94 y i 495346048 ont? . 
na mire 4 se wf 
- 
+ 
a Ac : P io 4 
4 a 
j 4 i J rh ne 


i Mr i 90m ereea 
Sid nraS oP) ™ pha 


i wert. | poe ihe met 6 ene ol 

r bb Sceeaaaee 0 uyeey 
4 | ivy || we > ie a at 
rt sbdom ong 


rk) | % 
| 
| 
La 5 
i 
; 
‘ "& 
‘a © 
4 ( 9 7 oP 
A 


“LISS a NF 
. , 
(SAy ! 
j a 


bo rioanys 
7 


ne 


aa 
oud 


ae 
an Ga 


sp 


A2 


is a measure of the Linear dependence of the two components 
xX. and Xj atwiagen, 

Characteristics of this matrix for 7 = ity are 
discussed in detail in the next section, together with a 
section on its frequency domain equivalent, the cross 
spectral matrix. A short discussion on the correlation 


matrix and cross spectral matrix with r 7 iP is included 


in Appendix Al.4. 


Al.2 The correlation matrix with r, = (Pas 6 applications to 


vectors which are restricted to a plane 


First consider the correlation matrix with ry = Fo: 
A symmetry property of the matrix can be readily shown: 


64 g(t) = <xz(t)xg(ter)> = <x, (tte) xg(t!)> 


<xy (tt) xq & fre)> bmp atin ae Oma eee 


Sle Mote (A3) 


Equation (A3) shows that 9; 5 (0) is Hermitean. 


Discussion here is limited to two very simple 
vector time series, totally random (or unpolarized) vector 
series and vector series restricted to a plane. Other types 
of vector time series lead to correlation matrices with 


very complicated characteristics [Roman, 1959]. In unpolar- 
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ized disturbances the end point of the disturbance vector 
moves quite irregularly. The vector shows no preferential 
directional properties and no predictable temporal proper- 
ties. Thus the spatially orthogonal vector components must 
always be uncorrelated Borey Saaaleeie = o|, and in 
order that temporal variations be purely random, <x. (t), 


* 
x; (t-t)> = 0. The further requirement that the terms 


TF#O 


5 
<x; (t), x (t-1)> be invariant under a coordinate 


ead 
transToymetion indicates that thescorre lation matrix UP of 
an unpolarized vector series is a scalar. More precisely 
U = Adi a> where A is a constant proportional to the power 
of the unpolarized disturbance. 

In a plane polarized vector time series the vector, 
when measured at one point in space, is restricted to a 
given plane. (Note the departure from the usual definition 


in optics.) Thus there must exist a lag-independent rota- 


Ghoneimacnixo | Rk Such that 
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a ape teat (A4) 


where #(t) represents the correlation matrix at lag t. 
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This rotation matrix can be obtained by diagonalizing 
6(0) [Cummings et al., 1969}. 

Consider now the correlation matrix (0) of a 
general vector time series. Since this matrix is Hermitean 
there must be a real transformation R which diagonalizes 
the matrix. The matrix in the new coordinate system, the 


principal axis system, is given by the transform 


AON 
RLo(0) IR! = | 0a, 0 (A5) 
OW Ne 


where we have chosen Ag < ho < hy: These eigenvalues are 
real and positive for any set of real time Series. Tne 
trace of the matrix is invariant under coordinate trans- 
formations and represents the total intensity, or total 
power, of the vector disturbance. The trace elements are 
the intensities in the individuahedirections.. Ihe matrix 


is conveniently expanded in the form 
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If we know dpriori that the polarized part of the vector 
time series is confined to a plane and assume that the 
polarized and unpolarized components of the pine Seves 
are uncorrelated, then the above expansion has a simple 
interpretation. (A) is the matrix of the unpolarized 
vector time series while (B) and-(C) are the two spatial 
components of the plane polarized vector time series. The 
eigenvectors corresponding to hy and ho define the 
plane of polarization. If we require (1 ,-A3) and 
(A4y-A3) to be positive then this expansion is unique. This 
is a reasonable requirement since, for a real time series, 
these two values are the two orthogonal intensities of tne 
plane polarized vector series. 

In general the polarized vector series will trace 
some complicated 3-dimensional shape. One helpful (though 
not conclusive) method to determine whether the postulate 


Ofpiane polarizacionems  COVVect SecLondevermiune 


be o Rie ae (A7) 


Lf o; 3(t)? 0 (where A, = 44(0)) for any t #0 then 
the series are not plane polarized and the above inter- 


pretation of expansion (A6) is incorrect. 
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A useful parameter in determining the degree of 


plane polarization is given by the expression 


R, = ——>. + (A8) 


Here Ry is the ratio of the apparent plane polarized 
Sugnal intensity to the total signalmintens ity sal Ry 


equals unity then the signal is completely plane polarized. 


Kies she cross=spectrad= matrax: applications to plane 
polarized vector fields 
The information contained in the correlation 
matrix can also be expressed in terms of a cross spectral 


matrix. We use the Wiener-Khinchin relations 


(ee) 


ste, aI bea (0) CoamMME dt (A9) 
o;5(t) J Sauna df (A10) 
GOrnind 4 


o; ;(0) | S,,(f)df (Al1) 


The terms in (All) may be considered to be the total auto 


and cross powers of the vector disturbance. Suppose now 
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that the disturbance has a spectral width Af and mean 
frequency ee Then we can estimate the correlation functions 


o; (0) from the equation. 


iN Ar 
Ue ala) ‘| S; ;(g)dg (Al12) 
f-Af/2 
In particular we have 


We could alternatively use the Fourier-Stieltjes transform 
to formulate this function [Wiener,1930]. It can be shown 
[fenkins and Watts, 1968,p467] that the matrix formed from 
the terms (A1l2) is a positive, semi-definite Hermitean 
matrix. As "Af >= 0, 44) ~ Sue where $4 5'7) are 
the auto and cross power spectral densities. 

To find the apparent plane of polarization (if 
one exists) we diagonalize the real partyof the cross 
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In the expression (A14) R; J, and. J , are in general 


° ” i] 
functions of f. For convenience we have chosen J33 < 


Joo < Jay: ihe transform of thertovalematrixne1s sthen 


J4] Unione © ue 
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Rk? Gaia: lok al Danan (A15) 
eee eo J 33 


Once again if J33 = Jn = J13 =05 the svector 


Serves: 1Seatruly apiane polanized(i.ea restnictedmstomdsplane). 
i} i] 

Lae OnaLy J13 and J53 are zero and if the non plane polar- 

ized vector series iS approximately random or unpolarized, 


then the cross spectral matrix of the plane polarized 


vector series is given by 


(044-933) he : 
cee aad 0 (A16) 
0 0 0 


TiseMali Xe Sed OUNGSDY cdi exXPanSmsOnmeSuilimnd imc Oma tain 
equation (A6). In analogy with (A8) we can also define 
a polarization ratio for the terms in the cross spectral 


matrix. This ratio is given by 
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R, @1y = 933) + ze - 933) (A17) 
Jas va von 55 


Following the work of Born and Wolf [1959] and 
Pow Lense tal. 5 | 1S07s|) 11 san 0SS ib emvons MOwmtha teste 
plane polarized part of the vector time series is composed 
of a vector which moves randomly in the plane and a vector 
series in which the two orthogonal components are completely 


correlated. We can accordingly expand the matrix (A116) in 


the form 
Ce ee ren Lv Tie eo 
= J $ (A18) 
ee (Gan Ia) Org -iPin Poo 


where for simplicity we have dropped the zero terms of 

the matrix. U is the matrix of the random plane polarized 
component, and P the matrix of the completely correlated 
vector series with det P = 0. Consideration shows that 

U is simply an eigenvalue of the matrix J, although the 
Unitary transTormatvom is) not, "In general... amted msrota tion 
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We choose the smallest value of U since the other 


gives negative values for Pay and Poo: 


Another useful parameter, the ratio of the 
completely coherent signal intensity to the plane 


polarized signal intensity, is given by the expression 


é F 1] 
(2 : 34 ‘ (05, . 3) 


Note that this expression is meaningful only if Ji = 0 


or if the non plane polarized vector series is truly 
random. 

When Af/f -<< 1 the vector series has sinu- 
soidal spatial components and the matrix P is simply 
thescross spectral matrix of an el liptically polarized 
vector time series, with the coordinate system oriented 
along the Major andemi nor aXxesu Ot mulicme ll ipSenesshile 


rotation eigenvectors of R_ corresponding to Jaq and 
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respectively. Ihe sense of rotation on sthe vectors right 


J are along the major and minor axes of the ellipse, 


or left handed, is given by the sign of Jiro or Pao: 
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All 
The total intensity of the coherent Signal is given by 
the equation 


(A22) 


Often a vector time series will have a band- 
width a Such that = 2 ei hd Ate oe Srl. Carn 
this case we might say that we are looking at the spectral 
components of the series and the associated frequency 


dependent ellipses of polarization. 
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where 77 (%,,77 0) is the power in the ith component 
at position re These functions have the simple 
characteristic 


Jog s(t) | | (A24) 


thesparametern  p is normal lyeanconp li catedurunctionmcr 
the lag t and is, except in some very simple cases, 
difficult to interpret. A detailed study of the cor- 
relation matrix has been given by Roman and Wolf [1960a, 
1960b] and Roman [1961]. 

The elements of the cross spectral matrix are 


defined by the relation 


ftat/2 


> > > > 
J5()-%o>f) 4h S;5(1).%.9) dg (A25) 
eee 


This equation is simply a generalization of equation 
(A12) with +, # r,- All the terms, including the 
diagonal elements Jay Joos J33 are in general complex. 


More precisely 
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twomseries® Itlis iconvenientetondersinessin analogy 


with equation (A23), a coherency estimator given by 
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The coherency kK can be considered as a measure of 
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the correlation of the two series at each frequency f. 


The phase spectrum estimator 
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is another useful parameter. Fag shows to what degree 
the frequency) components in one’ series lead or lag those 


at the same trequency in tie oOther=series: 


wiaen ver 
(rea Lay ligFi*ss y+ (%, aie mi US My of 
Oe Pin Viggt ararad cs 7 1 . 


nD _ 


a 

10 - Syn dist © eaober sh ee ed Mbp fe’: meh anes 5a 
J. omRdPS a Vea) fp ee eee er” Ome, ‘% vatsinyad, a fe 
ms ‘o - 


int emeese Suntesee veong aT mA, -. 


* CP, 95.0 7 O- : md), : wi 
fb ; — mane ~~ ; : 
ORK) TERE RS neg 4K) » Pgh NETS © 


’ = 7 bal ve a we © ae ron es ~ er ne wee 
S$yG00 244wW at (node 53" «12 7ONH OG i taed prrery wii : 


deo) gels hee! ostrab sag wh Ba neneanes Voraunent 243 
_isivee verte s@t_cb yeesupan? emake ‘ody Te 


Al4 


APPENDIX A2 ESTIMATION OF THE COMPONENTS OF THE CROSS 
SPECTRAL MATRIX; THE POWER SPECTRA AND 
CROSS SPECTRA | 


In order to obtain accurate estimates of the 


powers and cross powers J (equations (A12),(A25)) 


ds 
we need a spectral window Q which is approximately 


rectangular and with spectral width Af chosen such 


f tat/2 
QCF)*S(F) ig = sf he (A29) 
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where tne operator * denotes convolution. In our analysis 
the values Jig were calculated by smoothing the periodo- 
gram estimates [Jones,1965] which gives (as will be seen 
later) an approximately rectangular spectral window. 

The periodogram estimates for the power spectral 
density are obtained from JA(n)|° where A(n) is the 


discrete Fourier transform of the time series Xy(t). 
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transform algorithm (FFT) [Gentleman and Sande,1966] 

is used in the computation of A(n), Xv normally 

consists of the original time series X (m Werbaes jie 

with zeroes added to give N_ points (N being an integral 
power of 2). The periodogram estimates corrected for 


the added zeroes are then given by N |A(n) |¢ 
= : 


It is possible to show [Blackman and Tukey, 1958] 


that 
ne JA(n)|* = a Ja, (nat) |° (A31) 
where 
ay (f) = >, a(ftkat™ |) (A32) 
keno 
and 
a(f) -f aca) cou oe (A33) 


If there is no aliasing, which we now assume,then 


ay (f) = BP) Ome fia gis 


We obtain the time series Xn from the original 


signal x(t) by multiplying it by a data window eu) S 


hoa MGs) Bt) aie) (A34) 
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| B(t):x(t) e7°™' © at | “= p(F) (A35) 
WAG 


where A(n) is the function evaluated at nAf and es 
=mAt, If the signal x(t) is stationary, we can prove 


[Blackman and Tukey, 1958] that 
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S, the true power spectral density, is given by the 


relation 


Speen ay) ea) ee a (A37) 
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The spectral window W(f) is given by 
©0 Z 
W(f) = 7 B(t) e oT tt at (A38) 


To reduce the variance of the individual power 


2 
estimates jaca) |®, we determine the smoothed spectral 
estimates J from the equation 
k 
2 2 9 
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p=-k, 


The expression (A41) gives the estimates of the total 
powers UG in equation (Al2). The overall spectral 


window is k 


Q(f) = 7 a W( f-pAf) (A42) 


p=-k, 


We will show later that this window has an approximately 
rectangular shape if B(t) is properly chosen. 
The unsmoothed cross spectral estimates of two 


timesseries X(t) and Y(t) arewoiven by 
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neon) (A43) 


where A(n) and C(n) are the discrete Fourier trans- 
forms of X(t) and Y(t) respectively. To obtain the 
smoothed cross powers to be used in equations (A12) and 


(A25) we average the real and imaginary parts of the 
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Spectrum separately. The spectral window in each case 
is identical to that given in (A42). 
In our analysis we have used the time domain 


window B(t) shown in Figure Al. 


-Q-b =5 ° b atb 
Fig. Al. The data window B(t). 
B(t) = 5 (cos (S"(t+b)) + 1) yee 
27 
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= ] ala ee he xe Gy (A44) 
= 0 otherwise 


The transform of this window is 
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Some simple limiting cases are 
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T 
[Jenkins and Watts, p 244] 


lim 7 = 2b _sinwb 


a>0? ab (A47) 


[Jenkins and Watts,p 244] 


To find the shape of the spectral window Q(f) 


we must substitute W(f) =. |T(#)|* 


in equation (A42). 
Some typical spectral windows are shown in Figures A2, 


A3, A4. 
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‘In our analysis of the magnetic data, the window 


B(ty™ was set with a/b = O22. and the summation limits 


3N 


were k y= Ko= k= ae where k is the smallest integer 
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estimates (see equation (A39)). For convenience, the 


corresponding spectral window is denoted by 
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APPENDIX A3 VARIANCE OF THE SPECTRAL ESTIMATORS: 
CONFIDENCE LIMITS 


If a time series has an approximately norma | 
distribution of amplitudes, then the real and imaginary 
parts of the discrete Fourier transform of this series 


also havesnormals distributions eaeLneOLiereworasei] ft 


IMGT Sto eR) eae) (A48) 


where A(n) is the discrete transform, then a(n) and 
b(n) have approximately normal distributions [Jenkins 


and Watts, 1968]. The periodogram values 


JA(n)|* = a@(n) + b@(n) (A49) 


have x5 distributions and the smoothed periodograms 
; Z : : F , 
(equation (A39)) have X2 (2k+1) distributions with 
accordingly decreased variance. 
The variance V of the spectral window Q(f) 


(equation (A42)) can be determined from the relation 
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v= 8 nf a°(a)dg (A50) 


[Jenkins and Watts, p241] 


where S(f) is the power spectral density at frequency 


f, T the length of the sample, and Q_ the normalized 
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spectral window ince) Q(g)dg = 1). If we assume 
that the window Q is approximately rectangular with 


Width {2keW)At = (2ki1) (NAt\s othenathe wvartancelts 


yeanleSalt) (A51) 
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Ihe calculation of thesdistribution, functions 
and associated variances of the parameters obtained from 
the cross spectral matrix (see for example equations (Al7), 
URZOEM AZ.) aa Doles Ned, ee Aco eniSmO Gert damental 
importance in assigning confidence intervals to these 
parameters, especially since the expectation and variance 
of the parameters depend on the spectral window Q(f). 
It is particularly useful to know the distribution functions 
of the polarization ratios Ry and Ro and the coherency 


_K in order to eliminate those estimates which might 
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Gorres pond eto ta prandomenoise sourceseninewcalculaliongor 
these distribution functions is not an easy task since they 
are all non-linear functions of random variables. 

An approximation used to assign confidence limits 
to the phase and coherency spectrum is discussed by Jenkins 
and Watts. They have shown, using a Taylor series expan- 
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(A52) 


where I -f a*(F) dF, T the record length, and Kay the 


expected coherency. The transformed variable Y = tanh! 


Kay 


has an approximately normal distribution with a constant 
variance I/2T. The confidence limits to be used for the 
rejection of unwanted random noise can then be estimated 
from the variable Yos transforming back to the original 
scale to find Kig: WeStind. Us iINGmthismmecnodee tie tard 

95 per cent confidence limit for the coherency of random 


noise is 


Og< Kee. 


ij £0.48 (A53) 


where we have chosen the spectral window Q(f,k=3). 

The distribution functions for R, (equation 
CARZ a eand Ro (equation (A20)) were generated empirically 
by producing three independent, random series on a computer. 
PIcuresASmiss denis COgral sof athe distri DUtOnMOr R, obtained 
using the spectral window Q(f,a/b=0.2,k=3,N/m=1), and 
Figure A6 is a histogram for Ro obtained using the same 
spectral window but with the assumption R, = 1.0 (only 
two non-zero time series were used). Inspection of these 
ficures shows) that in Doth casesmlomapenscent Ofe thanestimates 
have values less than 0.75. Both the ratios R,, Ro=0.75 


can then be used as 95 per cent confidence limits for the 


rejection of random signals. 
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Fig. A5. Histogram of the values of R, which were 
calculated using the spectral window 
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Fig. A6. Histogram of the values of Ro which were 
calculated using the spectral window 
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APPENDIX A4 PRACTICAL CONSIDERATIONS IN ESTIMATING 
THESSPECTRALGS PARAMETERS 


A4.1 Selection of the data windows 


In choosing the length and shape of the data 
window B(t) (equation (A44)), a number of conditions must 
be met. First the lowest frequency and the spectral sample 
Spacing are set by the length mAt’ of the window. More 
precisely the spacing of the estimates is Af = (mat)~! 
(equation (A30)). Adding zeros to the time series is useful 
if the data have little noise since the frequency spacing 
is then (Nat)7! which allows for samples at lower fre- 
quencies. This does not increase the resolution of the 
spectral window which actually depends on the shape of B(t). 
All that the addition of zeros does is to give a finer spac- 
ing of the spectral estimates and in some cases the finer 
Spacing may actually be helpful in selecting very sharp 
spectral peaks. If the data are noisy, which is the case 
for most micropulsation data, then the spectral estimates 
must be smoothed in the manner shown in equation (A39) to 
reduce the variance. Since the variance of a spectral 
window is proportional to ‘ (equation (A51)), the estim- 


ates obtained from series with added zeros must be smoothed 


a proportionately greater amount to achieve a variance 
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equivalent to estimates obtained from series with no added 
zeros. The end results, whether zeros are added or not, 
are spectral windows with the same resolution. One final 
condition that must be considered is that all discrete 
Fourier transforms obtained by the FFT must have a number 
of points which is an integral power of two. This con- 
straint can place severe limitations on the length of the 
sample so that one normally adds zeros to make up the 


desired number of points. 


A4.2 Detrending the data 


It is sometimes apparent from inspection of spectra 
that mucn of tne power is Contained im &@ Few Very Harrow 
bands, especially at low frequencies. Since spectral 
windows have some leakage (see for example Figure A3), these 
peaks cause errors in spectral estimates at lower powers. 
Because of this, it is usually necessary to digitally 
filter the data to remove power at the frequency of these 
bands before spectral analysis 1S» Carnieds outs ins magneric 
data recorded with fluxgate magnetometers, low frequency 
trends can have amplitudes hundreds of times larger than 
an associated micropulsation event. Consequently it is 


advantageous to use filters which supress the low frequency 


components of the magnetic signals. 
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A4.3 Characteristics of the data windows and detrending 
Tall cerns 


In the analysis of the magnetic data, two differ- 
ent data windows were used. The first was 60 min long 
and was used in conjunction with a 1-20 mHz bandpass, 
zero-phase shift, digital filter [AZpastan, 1968]. The 
second was 30 min long and was used in conjunction with a 
2.5-100 mHz digital filter. The frequency and impulse 
response functions of these filters are shown in Figures 
A7, A8, and A9. To eliminate the transient response of 
the filters the first.sand last J0°minwot fisiceredsdata 
were rejected. Both data windows had a cosine taper with 
a/b = 0.2 (equation (A44)). This cosine taper was used 
to eliminate any remaining traces of the transient response 
of the filters and to improve the shape of the smoothed 
spectral window. In the analysis of the 1969 data only the 
60 min data window was used and the data were sampled at 
every 28th point. In the analysis of the 19/70 data using 
the 60 min window, every 16th data point was taken, and 
in using the 30 min window, every 8th data point. This 
gives a total of 117 data points for each window with the 
ratio N/m = 128/117. Each smoothed spectral window Q(f) 
had k=3 (equation (A39)) with spectral window widths of 
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1.78 mHz and 3.56 mHz respectively. With filtering and 
spectral smoothing, the 60 min data window provides a 
useful frequency range of 1.0-15.3 mHz,while the 30 min 
data window provides a range of 2.5-30.5 mHz. These 
ranges correspond roughly to the Pc5 and Pe 4 micropul- 


Sation bands. 


A4.4 Selection of spectral peaks 


The selection of spectral peaks: 1s euditticult 
problem since the spectral window Q(f, k=3) has very 
poor resolution at low frequencies. If we define the 
resolution of the spectral window by the ratio Af/f where 
Af is the width of the snectral window, then we find, 
for example, that in using the 60 min data window af =A lao 
at a frequency f = ImHz. A solution to this problem is to 
estimate the positions of the spectral peaks from the un- 
smoothed power (equation (A30)). A better method is to 
smooth the power estimates (equation (A30)) using a window 
Q,(F) with width Av(f) (Av being a function of f) such 
that 


AU = (A54) 


where R is the best resolution desired. For the inter- 
pretation of the polarization parameters it does not matter 


that the window Q,(f) changes its spectral width and 
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Variance across the spectrum, since it is used only to 
locate the spectral peaks. We can then go back to the 
Original spectral data, which was smoothed Tite GZ (ess 
to determine those estimates for which the window Q(f, k=3) 


covered only one spectral peak. 
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APPENDIX A5 THE EFFECTS OF DIGITIZATION ON THE SPECTRUM 


The digitization of data adds errors to the data 
field which depend on the discrete time sample interval 
At and the amplitude sample interval Aq (henceforth 
referred to as the quantizing interval). Under simple 
constraints it can be shown that the r.m.s. noise error 
arising from digitizing is given by (Aq)*/12 [Bendat and 
Prersoe. 91960 |; pine tactwit i1sepossib lestompnovemcnac, 
independent of the data signal used, there is an upper 
Dound= to, the: digitizing, erroremt stiemd ldic12inosmolsens.s 
assumed to be a white Gaussian noise with total power 
(Aq)*/12 and bandwidth (2At)~! EMazo, 1968]. In practice 
the digitizing error often has the same effect as an added 
noise source and the spectrum of the digitizing error is 
independent of the signal over a wide range of signal 
amplitudes. Even if the signal is so small that its 
amplitude covers only several quantizing intervals, there 
is usually enough background instrumental and magnetic 
noise to determine the amount of digitizing noise and hide 
its relation to the signal. 

More precisely, Bennett [1948] has shown that if 
the digitized signal is a random noise function with a 


frequency spread 0 to lees then the power spectral density 


a 


eh bes ae ab 


eadtaate 4uT dh’ wath’ We re ae 


oe 


eob SAY oP e@erne4e atts e76b Fe, worden? fb i 

[< oy 7 

‘nsetet afomes onit aeioeth aie na’ tnngisb’ to Yds 
= 


a) 
‘9. $@Pan 2227.4 nie 20 A069 om as 2003 
| ey iy f patsttigss wor? eh sie 
Py 
acu t - 7707 wt’ « { saer Some att 
i 4 i? nied o43 TA a 
hehe ote: yd” be 7 


Di d hy 
Te 70799 Oni ran 
ao bae sme Od 
mye Yo | bee s teva Jeeste e627. 40 tnabagk dole 


laws oe 4 Tom iy 3) need aabut Hg 
Af. rie: rs 164 2 oi Me 15 fo UP Yon 7 
4 el 12y¥ ayaa + 7 


otiadper See Catered Ant mnpdigdeed | Agvons r pads zt 

: ib fon wart me ¥ J arianete oF s2to 
ieedibe or of a0 otestan 

+t od inwode ean Labbe) ieepael d qo 909 * om 


a «a an ‘Ze baxts: 


n m7 tu.eotissut s eon nob ’ 
; . : 
i 8 ob aot nt 4 > ; 


y279 nen lerioacge Yoweoq, ane rod . 
= 4 6 - 
- 


a as 


A33 


Q(f) of the quantizing error, with no time sampling, is 


given by 
we 
Ue eae 3(aq)efe| 
PACE Ne Sy ARES EECEE Es gh SLCAU LE (A55) 
Ono f 212 Bria f° 
0 0 econ 6 
where 
= ooxin’ 
B(x) = > aaa (A56) 
n=1 p 


and 25 iS the 7.ms. SignalsValucem nla werinclude 
discrete time sampling, then the power spectral density 


function, with aliasing, becomes 


[ee] 


A(t) = SS 0(f+2nfy) 0<f<fy  (A57) 
n=-0 

where fy = ar ae the Nyquist or alias frequency. 

If the spectrum of the quantizing noise is aliased many 
times, A(f) becomes nearly uniform between the Nyquist 
frequencies. Since quantizing the signal amplitude 
effectively spreads the bandwidth we would expect consider- 
able aliasing of the quantizing noise. With this in mind, 
we find that the aliased power spectral density is approx- 


imated by the expression 
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In our experiment f 
spectral density of the quanitizing noise is 0.17y°Hz 
If an input signal is quasi-sinusoidal with a peak to peak 
amplitude of 2Ay, then the ratio of the signal power to 


the noise power is ones 
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APPENDIX A6 A BRIEF DESCRIPTION OF THE INSTRUMENTATION 


The basic components of the magnetic recording 
System are shown in Figure Al0. Each station "Ee equipped 
with a three-component magnetometer combined with a WWVB 
receiver, analogue-to-digital system, and digital tape 
recorder. Data were recorded on seven-track magnetic 
tapes at a rate of 3.13 samples/sec in 1969 and at a rate 
of 1.56 samples/sec in 1970. The records were timed by 
recording the WWVB radio signal directly on tape every 
lophrecs minoan 1969 ‘andeevery 7 wireocemine inethemio70 
field season. 

All magnetometers, except tnat at Leduc (where a 
Sharpe MFO-3 was used), were three-component fluxgate 
magnetometers designed by the Geomagnetic Division of the 
Department of Energy, Mines and Resources. The theory 
behind the fluxgate system has been discussed by Sersen 
[1957]. Trigg et at.[1970] give a complete outline of the 
type of magnetometer used in this experiment. The latter 
paper provides the general design information as well as 
details necessary to operate the magnetometer correctly. 

The component parts of the magnetometer are 
displayed in Figure All. In our experiment we mounted the 


sensing heads on aluminum poles set in concrete at the 
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Fig. All. The components of the fluxgate magnetometer. 
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bottom of rigid plastic containers. The containers were 
buried with their tops at least 10-20 cm underground. 
The outputs (X, Y, Z) from the amplifiers had an optimal 
sensitivity of 10mV/y and were linear over most of the 
amplifiers' +10 V range (Figure Al2). The gain and phase 
shift of the complete magnetometer are essentially flat 
to 0.1 Hz (Figure A13) which is well beyond the range 

of frequencies of interest in our experiment. 

To prevent aliasing in the spectral analysis of 
the data, all magnetic signals were passed through two- 
Stage low-pass filters. Even though the power in the natural 
magnetic spectrum falls off rapidly with frequency, there 
is often enough natural and instrumental noise at high 
frequencies to permit considerable aliasing. The frequency 
response characteristics of a typical filter are shown in 
FiguiresAl4a) lteeaddi tion, thesdliias ingen! itersaproviudesgaod 
protection against voltage overloads to the multiplexer 
Since the maximumerange Of —tN@ea il lensmiS = |i2avOlts. 

The analogue-to-digital system is based on the 
Redcor 720 A-D converter, with an added multiplexer bypass 
line for WWVB, an added WWVB-to-data relay, and a variable 
stepping rate system coupled to a crystal oscillator 


(Figure Al5). The crystal oscillator and its associated 
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Fig. Al3. Gain and phase response of the magnetometer. 
(after Trigg et et. [1970]) 
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divide-circuit control the sample rate of the system. 
Alterweacheset a7 /#data Dlockse( /anrescemin) thes stepping 
rate changes from 1.56/sec (corresponding to a sample 
interval of 1.92 sec/component) to 50.0/sec in order to 
permit an accurate recording of the WWVB time signal. 

Each recording of the WWVB signal lasts 123 sec. Since 
the maximum drift of the oscillators is + 0.005 per cent, 
the absolute timing should not be in error by more than 
two sec over the interval between WWVB recordings. (In 
1969 the interval between WWVB recordings was 16 hr 23 min, 
and the sample rate was 3.13/sec). The sequencer controls 
the sampling order in the multiplexer. The samples follow 
a cyclic order (HDZHD...) and all the sample intervals 

are of the same length. The 12 bit A-D converter uses 

the successive approximation technique with an offset 
binary digital coding format. The logic format board 
arranges each sample into two bytes, each byte being 6 
Ciugicalebl tS) Ond mine meaS tes TOlinalcal ua) ktm oma LODDoG 
Fomperni t this Ditepos iti One lOmDes Us edmccmanchanmelen ag 
(the flag is utilized to identify errors in the sequential 
sampling of the three input channels). Considering the 
full scale input for the A-D converter (+10 volts, cor- 


responding to +1000y) and the available 11 digital bits, 
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we find that the maximum range of the system is +1000y 
with an amplitude sample interval (or quantizing interval) 
oP Ml00072 =) 0.9767. 

The output from the analogue-to-digital system 
is recorded on a seven-track Peripheral Equipment incre- 
mental write tape deck. Each sample is two bytes long 
(2 bytes x 6 bits) and each data block is 6144 samples 
long (12288 in 1969). This block length is equivalent 
LOLOD min sUssec. of Magnet 1 cCedatasandmico msec or stne 


WWVB time signal. 
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APPENDIX AZ EARTH INDUCTION EFFECTS 


In any attempt to relate ground-based magnetic 
fields to their sources in the ionosphere and magneto- 
Sphere we must assess the effects of induced earth currents 
and eliminate such effects from the data. In the discussion 
here we outline briefly the induction effects of uniform 
and simple, non-uniform earth models on low frequency 
electromagnetic disturbances. 

Price [1950, 1962] has shown that the magnetic 
field at the surface of a uniform conductor (Figure A116) 


is given by the equation 


H= - |(a+8) 2 a) sve (-AeB)] ef 2"F (A59) 
where the inducing field is 
Micah [ ead, vp] e124 (A60) 
P is the magnetic potential 
B= 2 A (AG1) 
and 
e@ = y* + i2mifo (A62) 
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The reciprocal of the parameter v is a measure of the 
NOrEZontdales cCaleeo thes OUrGemnlel Gm \e=s 217 vieiate 1s 
the source frequency, wp the magnetic permitivity, and 


o the conductivity of the ground 


Fig. Al6é. The coordinate system for the uniform earth model. 


The ratio . = ed gives a good measure of the relative 
induction effect. 

Inspection of equation (A59) shows that uniform 
induction has no effect on the horizontal projection of the 
polarization ellipse. The case is different, however for 
the x-z and y-z planes. As induction increases, the 2 
component decreases and changes phase and the major axis 
of the ellipse tilts toward the horizontal plane. For 


complete induction (o + ~) the z component is zero and 


the two horizontal components increase by a factor of two. 
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TablesAl illustrates the relativesindcuctiom effects tor 


ValouseValues of 9G 2nofuav = {10 | Joti. 


TABLE AL. Relative Inductionsbtfects 


Evidently induction becomes noticeable only when 8B > 1. 


If we take the range of values 


Ga aoc 


10 
Weed coun ak 3(107°)Hz (A63) 
‘War < "gr < i awe fa 


then the range of 8 for the parameters influencing our 


experiment is 


NACE Re RS ELS) (A64) 


These limits show that induction will have an appreciable 


effect except at low frequencies (< 10 mHz) in resistive 


regions (o < 10me) 
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The evaluation of induction effects near horizontal 
inhomogenieties presents a more difficult problem. In 
general the magnetic and electric fields are determined 
by numerical, analysis. (see, GudoJones,, 1970s suzft, 197.14. 

It can be easily shown [see e.g. Jones, 1970] that near a 
Venta calecontact’ ot contrasting resist weities Figures nds) 
the z component of the magnetic field is markedly en- 
hanced when there is a component perpendicular to the 


Surike. (1... 11 Y. GCIrectron: 


Fig. Al7. The coordinate system for the model considered 


by Jones [1970]. 


Schmucker [1964] noticed this effect in an analysis of 
geomagnetic data near the coast of California. The relative 
phases of the y component and z component are seldom 
greater than 20° in this model and thus the polarization 


in the y-z plane is nearly linear in all regions. The 
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X Component remains constant across the strike because 


of the two dimensional character of the problem (i.e. 


Rankin and Reddy [1971] have considered a rather 


more complex earth model (Figure A18). 


Fig. Al8. The coordinate system for the model considered 


by Rankin and Reddy [1971]. 


Figures Al19-A21 show the effects of the discontinuities 
in conductivity in this model on the horizontal magnetic 
field components. The application of their model is 
limited to sources of infinite extent but the model does 
GiVemcaecuCOdinclcauiun Ure cine sIndUCUlOne@ehLeGtl Old) > — 
continuities similar to those found in coastal regions. 
The greatest induced changes in the sense of polarization 


occur when the source field is linearly polarized at an 
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Fig. Al9. The relative amplitudes and phases 
of the y component of the magnetic field (Hy) in 


the model considered by Rankin and Reddy. ty is 
the y component of the magnetic field at large 
distances from the vertical contact. (Note the 
change in scale on the abscissa; after Rankin and 


Reddy [1971]). 
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angle of 45° to the sunhiea (0Phente A20)2 These lip tici ty 
changes from ‘0 on the resistive side (land) to ‘0.1 
on the conducting side (ocean). If the source is cir- 
cularly polarized the apparent sense of polarization 
remains unchanged (Figure A21) although the magnitude of 
the ellipticity varies from 1.0 to 0.5. “Note that the 
MajOreaxiS of then ela ipse 15 Ddaleclasco. tess tra kKeeon 
the resistive side and perpendicular to the strike on the 


conducting side. 
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